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Abstract 

A review of selected topics in Hirota's bilinear difference equation (HBDE) is given. This famous 
3-dimensional difference equation is known to provide a canonical integrable discretization for most 
important types of soliton equations. Similarly to the continuous theory, HBDE is a member of 
an infinite hierarchy. The central point of our exposition is a discrete version of the zero curvature 
condition explicitly written in the form of discrete Zakharov-Shabat equations for Af-operators 
realized as difference or pseudo-difference operators. A unified approach to various types of M- 
operators and zero curvature representations is suggested. Different reductions of HBDE to 2- 
dimensional equations are considered. Among them discrete counterparts of the KdV, sine-Gordon, 
Toda chain, relativistic Toda chain and other typical examples are discussed in detail. 
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1 Introduction 



In 1981 R.Hirota published the paper H] which summarized his earlier studies of discretizing nonlinear 
integrable equations The main result is a compact bilinear equation which can be viewed as an 

integrable discrete analogue of the 2-dimensional Toda lattice. In Hirota's original notation it has the 
form 

[Zi exp( J D 1 ) + Z 2 exp(£ 2 ) + Z 3 exp(£> 3 )]r • r = , (1.1) 

where Zi are arbitrary constants, r = t(xi, x 2 , x 3 ), Di = D Xi and Hirota's D-operator is defined for a 
linear differential operator F(d x ) by 

F(D x )f(x)-g(x)=F(d y )f(x + y)g(x-y) . (1.2) 

y=o 



In the more explicit notation eq. ( |l . 1| ) looks as follows: 

Z\t{x\ + 1,x 2 ,x 3 )t(xi - l,x 2 , x 3 ) 
+ Z 2 t(x 1 ,x 2 + 1,x 3 )t(x 1 ,x 2 -1,x 3 ) 

+ Z 3 t(x 1} x 2 ,x 3 + 1)t(xi,x 2 ,x 3 - 1) = 0. (1.3) 

This equation is often called Hirota's bilinear difference equation (HBDE). Its simplicity is surprizing and 
delusive at the same time: each detail is controlled by intcgrability and hides meaningful mathematical 
structures whereas some even simpler looking equations turn out to be untractable by analytical methods. 

One of the most impressive outcomes of Hirota's works is that HBDE is shown to unify many if not 
all known soliton equations. More precisely, it contains them in an encoded form. Performing a scaling 
continuum limit for appropriate combinations of parameters and variables, one is able to obtain the 
Korteweg - de Vries (KdV) equation, Kadomtsev-Petviashvili (KP) equation, modified KdV (MKdV) 
and modified KP (MKP) equations, two-dimensional Toda lattice (2DTL) equation, sine-Gordon (SG) 
equation, Benjamin-Ono equation, etc. Their discrete analogues are pro duced from HBDE by choosing 
suitable dependent and independent variables. Furthermore, eq. (fTl]) was shown to possess soliton 
solutions and Backlund transformations for generic values of parameters. These facts suggest to consider 
HBDE as a fundamental classical soliton equation, from which the typical examples can be obtained as 
particular cases. 

Recently, bilinear equations of this form emerged , § in the context of quantum integrable systems 
as the model-independent functional relations ||, Jl0| for eigenvalues of quantum transfer matrices. This 
was our motivation for revisiting classical nonlinear difference equations. 

These notes aim at reviewing selected topics in HBDE and further clarifying basic elements of the 
theory. In our exposition, we deal solely with equations themselves saying almost no word about their 
solutions^- Likewise their continuous counterparts, completely discretized nonlinear integrable equations 
are known to possess soliton and finite-gap solutions. However, a systematic treatment of these and other 
particular classes of solutions could be a separate enterprise which requires much more space. We shall 
confine ourselves to elaborating discrete versions of commutation representations and auxiliary linear 
problems on a formal algebraic level. At the same time some important elements of our approach are 
essentially motivated by the finite-gap theory. 

The difference soliton equations are intimately connected with the differential ones. We already 
mentioned that the latter are obtained from the former by a scaling limit. Better to say, HBDE was just 
designed to enjoy this property. The fact that such an equation does exist is by no means trivial. A link 
in the opposite direction was established by T.Miwa jll} who noticed that discrete Hirota's equations 
can be obtained from the continuous KP hierarchy by choosing the time flows to be certain infinite 
combinations of the standard flows of the hierarchy. This idea was further developed in the papers 
jjJI, jl3| as a method to produce discrete soliton equations from continuous ones. The interrelation 
between discrete and continuous integrable hierarchies looks like a kind of Fourier duality: they provide 
complementary descriptions of the same object, namely, of the infinite dimensional grassmannian 



1 Just because of this we do not draw any distinction between discrete and difference equations. Usually, in the latter 
case it is implied that solutions are functions of a continuous variable with certain analytical properties. 
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In this survey we do not give a systematic treatment of the connection between discrete and con- 
tinuous hierarchies. The problem to describe the limiting procedure that would be compatible with the 
whole hierarchy is technically involved. However, it is impossible not to refer to continuous hierarchies 
at all. We agree to a compromise restricting ourselves by a few typical examples. 

It is assumed that the reader is familiar with the basic notions of the continuous theory such as Lax 
and Zakharov-Shabat equations, zero curvature conditions in scalar and matrix forms, commuting flows, 
infinite hierarchies, r-function etc. 

Let us outline the contents of the paper. (More detailed descriptions are given in short introductions 
to each section.) 

Sect. 2 can be considered as a part of the Introduction. We tried to collect here different forms 
of 3-dimcnsional HBDE known in the literature. All of them arc equivalent. Simple transformations 
between them are listed. 

Sections 3, 4 and 5 form the main body of the paper. To figure out the key principles underlying the 
variety of integrable difference equations, we need a number of definitions and axioms. They are given in 
Sect. 3. All the notions explained in this section are essentially used in the sequel. In Sect. 4, the discrete 
version of the zero curvature representation is presented. Filling in some gaps in the existing literature, 
we give explicit forms of the M-operators (realized as difference operators) for discrete flows. Sect. 5 
is devoted to various types of associated auxiliary linear problems. They provide a "linearization" 
of the original nonlinear equation. The related notion of Backlund transformation is discussed. The 
Baker-Akhiezer functions are introduced as special formal solutions to the linear problems. 

Sections 6 and 7 are more technical and might be interesting mostly for experts. They can be skipped 
without loss of understanding. In Sect. 6 it is explained how to extend the M-operator approach to ar- 
bitrary discrete flows defined in Sect. 3. In the general case the M-operators contain negative powers 
of first order difference operators. In Sect. 7 we dwell upon hierarchies of bilinear discrete equations 
and suggest the notion of " higher" discrete flows with the corresponding zero curvature representation. 
We conjecture that all "higher" (AT-term) discrete Hirota's equations known in the literature are con- 
sequences of the 3-term ones. This assertion is proved for the first nontrivial example of the 4-term 
equations in 4 variables. 

Sect. 8 deals with 2-dimensional reductions of HBDE, the corresponding L-M pairs and auxiliary 
linear problems. The list of reductions includes discrete analogues of the KdV equation, the ID Toda 
chain, the AKNS system, the relativistic Toda chain, the sine-Gordon equation, the Liouville equation 
and some others. 

In the Appendix, we present main elements of a different approach to HBDE based on Miwa's 
transformation. This is the method for generating discrete soliton equations suggested in ref. fl2]| . We 
show how it works for the simplest examples and comment on the continuum limit which is in a sense 
an "inverse" Miwa transformation. 



2 Equivalent forms of the bilinear equation 

Hirota's difference equation exists in several forms. Historically, they emerged as integrable discretiza- 
tions of particular continuous hierarchies (e.g. KP, 2DTL). In this section we give a list of most popular 
forms of Hirota's difference equation and explicitly demonstrate that they are equivalent. However, it is 
useful to bear in mind all of them since one or another may be more convenient in a particular problem. 



a) Hirota's original form (1.3): 



Zxt{xx + 1)t(xi - 1) + Z 2 t(x 2 + 1)t(x 2 - 1) + Z 3 t(x 3 + 1)t{x 3 - 1) = (2.1) 

(here and below we often skip the variables that do not undergo shifts). Note that the 3 variables 
enter in a symmetric fashion and the equation is invariant under their permutations and a simultaneous 
permutation of ZiS. The equation is also invariant under changing the sign of any one of the variables 
and under the transformation 

t(xi,x 2 ,x 3 ) -» Xo(xi + x 2 + x 3 )xi(x 2 + x 3 - xi)x 2 (x! + x 3 - x 2 )\3(xi + x 2 - x 3 )t(xi,x 2i x 3 ) , (2.2) 
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where Xi are arbitrary functions. The transformation 

t(x-l,X2,x 3 ) 



r?-xy2 ri—xl/2 , , 

Z l Z 2 Z 3 t(x 1 ,x 2 ,x 3 ) 



converts eq. (1.3) into the canonical form, 

r(xi + 1)t(x! - 1) + t(x 2 + 1)t(.t 2 - 1) + t{x 3 + 1)t{x 3 - 1) = . 
This equation does not contain any free parameters. 

• a') "Gauge invariant" form: 

(l + Y( Xl ,x 2 ,x 3 + 1))(1 + Y(x 1 ,x 2 ,x 3 - 1)) 



(2.3) 
(2.4) 



where 



Y(xi,x 2 + l,x 3 )Y(xi,x 2 - l,x 3 ) 



Y(x 1 ,x 2 ,x 3 ) 



(1 + Y- 1 (x 1 + 1, x 2 , x 3 ))(l + Y- 1 (x 1 -l,x 2 , x 3 )) ' 

t{x 17 x 2 ,x 3 + 1)t(xi,x 2 ,x 3 - 1) 
t(xi + 1, z 2 , aj3)r(xi - 1, a; 2 , x 3 ) 



(2.5) 



(2.6) 



is a gauge invariant quantity: the "gauge" transformation (2.2) does not change it. This form is a 
discrete counterpart of nonlinear integrable equations written in terms of potentials and fields rather 
than r-functions. Some particular cases of this equation emerge naturally in thermodynamic Bethe 
ansatz fulfil. 



b) KP-like form: 



(z 2 — z 3 )t Pi+1,P2,P3 t Pi,P2+1 ' P3+1 



+ {Z 3 - Zl ) T Pl^ + l,P3 T Pl+l,P2,P3+l + 

+ ( Zl - z 2 )t Pi > P2 > P3+1 t Pi+1 > P2+1 ' P3 = , 



(2.7) 



where t Pi ' P2 ' P3 is a function of the 3 variables p^s and Zi are arbitrary constants. This equation is 
invariant under cyclic permutations of the variables and simultaneous permutations of Zi's. Changing 
signes of all the variables also leaves it invariant. Besides, the transformation ( |2.2| ) can be translated to 
eq. (2.7): if t Pi ' P2 ' P3 is a solution, then 

Xo (2pi + 2 P2 + 2p 3 ) Xl (2 Pl ) X2 (2p 2 )x 3 (2 P3 )T p i' P2 > P3 



(2.8) 



is a solution too. Again, the coefficients in (2.7) can be made equal to 1 by means of the transformation 

P2P3 

T Pi,P2, P3 } (2.9) 



T Pl,P2,P3 _ 




-z 3 y iP2 




-*\ 








\Z2 


- z 3 J 



bringing eq. (2.7) into its canonical form. 
• b') MKP-like form: 



(z Zl )(z 2 Zs)t%+ 1 *>*>t%% +1 *> +1 + 
+ (z z 2 )(z 3 - z 1 )i%*> +1 **t%#***> +1 + 

+ (Z Z 3 )( Zl - Z a )7%*>**+ 1 T%£*»+ 1 ** = . 



(2.10) 



Note that the combination of the arguments p\ + p 2 + p 3 — po is the same for all r-functions in this 
equation. In orther words, the hyperplane p\ +p 2 +p 3 —po = const is invariant. Therefore, this equation 
actually depends on three variables rather than four: say, p\,p 2 ,p 3 . Since sum of the coefficients in 
eq. (2.10) is zero, like in eq. ( [2.7| ), they differ by a reparametrization of Zj's only. 

• c) 2DTL-like form: 



(2.11) 
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where t„> is a function of the 3 variables and fi, v are arbitrary constants. The variables I, I are called 
light cone coordinates. Note that in this form the permutation symmetry is lost. However, an analogue 
of eq. m ) holds true: if t 1 / solves eq. then xo(2n + 2l)xi(2l)x2(2l)x3(^n - 21)t 1 / is a solution, 

too. The transformation 

—11 / x -n 2 /2 



J.J 



1 



J. I 



allows one to hide the coefficients in eq. (2.11): 



T l,l+X T l+i,l 



i,i 

n n 



1,1+1 i+i.i _ n 

T n+1 1 n-l — u ! 



(2.12) 



(2.13) 



which is refered to as its canonical form. 

For reader's convenience we present here the linear substitutions making the canonical forms of 
equations a), b), c) equivalent. 

a)<-> b): t(x 1 ,x 2 ,x 3 ) = T Pl ' P2 ' P3 , 



Pi 



\{-Xi +x 2 + x 3 ) , 



XI = P2 + P3 , X2 = Pi + P3 



Ux 1 -x 2 + x 3 ), p 3 = i(xi + x 2 - x 3 ) , 



b)^ c): tP^ 



T l,i 

1 n -i 



n=p 2 +P3i 
Pl=l, P2 



I = p x 

- n — I , 



a)« c): t(xi,x 2 ,x 3 ) 

n — x\ 



1.1 



I = h(-xi +x 2 + x 3 ) , 



Xl 



X2 



n + I — I . 



X 3 = pi + P2 , 

l=P2, 
P3 = L 

I = 7j(xi - x 2 + x 3 ) . 
X3 = I + I. 



(2.14) 
(2.15) 

(2.16) 
(2.17) 

(2.18) 
(2.19) 



Clearly, these linear substitutions are not unique. All other possibilities can be obtained from the 
given one by applying a transfo rmati o n of t he form (xi, X2, x 3 ) — * (iip^j, ±ip( 2 ), ±ip( 3 )), where P is 
a permutation. Using formulas ( |2.14| )-( 2.19 ) one can easily obtain gauge invariant forms of equations b) 
and c). 



3 Definitions: the nomenclature of flows 

Here we introduce a practical set of definitions and axioms which will help us to develop a system- 
atic viewpoint to the zoo of non-linear integrable equations and their commutation representations. 
This viewpoint is in fact more general than we need for HBDE itself. Differential as well as "mixed" 
differential-difference non-linear equations fit the scheme, too. Our approach is motivated by algebro- 
geometric solutions [^0| to soliton equations expressed through Riemann theta- functions. However, since 
the goal is to clarify formal algebraic structures, we never refer to the solutions explicitly. 



3.1 The variables and kinematical constraints 

The "unknown function" entering bilinear equations is always denoted by r. This function depends on 
infinite set of independent variables which are called flows or times. The last two words will be used as 
synonyms. For each particular equation only a finite number of the time variables take non-zero values. 

The flows are labeled by points of the complex plane C. Call points A € C labels. We make distinction 
between discrete and continuous flows^J. 

• Discrete flows: With each ordered pair of points A,/i £ C, A ^ /i, a discrete flow I = l\^ is 
associated. To put it differently, the flows are attached to vectors A/i, i.e., each discrete flow has 
two labels. 

2 These are not more than conventional names. In general the both time variables may take complex values. 
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• Continuous flows: With each point A G C an infinite sequence of times {t\, t%, £3, . . . }^ is associ- 
ated. All the variables tj have the common label A. 

In each particular equation which we are going to consider only a finite number of labels are involved. 
To axiomatize this situation, we say that for all but a finite number of labels A s C and for all but 
a finite number of ordered pairs of labels the corresponding variables are implied to be equal to zero. 
This condition makes the definition very close to the adelic ideology from the algebraic number theory. 
The definition in such an abstract form may seem to be overcomplicated and too general. However, this 
standpoint is useful since it provides an adequate formalization of the simple fact that the number of 
independent variables in equations of an integrable hierarchy can be arbitrary but finite. 

Sometimes it is convenient to say that those variables which are non-zero are swithed on while all 
others are swithed off. According to the above definition, the set of labels corresponding to the swithed 
on variables is always finite. 

Having this in mind, it is worthwhile to rephrase the definition making it a little bit more concrete^. 

Let {A a } , a € /, be a finite set of marked points in C. Here, I is just the finite set of labels 
corresponding to the variables that are swithed on. By l a p (a ^ f3) denote the discrete variable associated 

with X a Xp. By tj , j = 1,2,... denote the continuous times associated with X a . The r-function is a 
function of these variables: 

r = r({l af3 };{tf}). 

Let Q be the graph whose vertices are the marked points (labels) X a , a € I, and whose edges are 
vectors X a Xp. The edges have orientation that is indicated by an arrow looking from a to /3. This graph 
will be refered to as the graph of flows. It encodes the kinematic structure of the equation. 

We stress that the only essential elements of the graph are vertices and their ordered pairs. All other 
graphical elements are introduced for convenience of the visualization. In particular, the vectors may 
intersect on the complex plane but the intersection points should be considered as not belonging to the 
graph. It is also worth emphasizing that the vectors are just convenient names of flows. They should 
not be mixed with "directions" of the flows in any sense of this word. 

The introduced variables are not independent. There are certain " kinematical" constraints imposed 
on them. 

The first group of constraints involves discrete variables only. The constraints arise when the graph 
of flows G has cycles. It is enough to fix the constraints for the following two cases: 

i) The 2-cycle: 

T(l a @ + 1, ha + l) = T\l a l3, l/3a) ■ (3-1) 



I at 



c 



3-v 



l/3a 



Informally, this means that lp a is identified with —l a p- 
ii) The 3-cycle: 



'(lap + 1) + 1) If a + l) — T\l a p, Ipy, Irya). (3.2) 



3 For each concrete example this unified notation is still not very convenient to work with. In the technical part of the 
paper this will be changed and simplified. However, for the sake of clarity and definiteness, it is better to introduce general 
notions and definitions using the unified notation. 
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The corresponding rules for longer cycles follow from these two. According to these rules, one can 
subsequently remove all the cycles and reduce the graph to a tree. The tree graphs correspond to 
kinematically independent flows. Formally, it is sufficient to consider graphs without cycles. However, 
introducing cycles sometimes makes the set of variables more symmetric though non-minimal. 



The last constraint describes the interrelation between a discrete flow X a Xp and the "adjacent" 
continuous ones (i.e., corresponding to the endpoints X a and A^): 

iii) Miwa's rule 

r(l af3 + l; tW,t<«)=T(l afl ;tW-[A /J -A ] ) t^) I (3.3) 
r(/ a/3 -l; t^\t^)= T {l a0 ;t^,t^ -[X a -X ]), (3.4) 

t {a) I a t {(3) 

A„ • ► • As 



Here r(t) = r({tj}) and the short-hand notation 

f(t± [z]) = f(h± z, t 2 ± ±z 2 , i 3 ± ±z 3 , . . . ) (3.5) 



is used for a function / of the infinite sequence of variables t = {ti, t 2 , ■ ■ ■ }. The second relation (3.4) 
follows from the rule (i) and the first one. 



The relations (3.3), (3.4) should be understood as formal rules which allow one to translate the 
infinite sequence of continuous time shifts into the shift of a single discrete variable and vice versa. We 
do not care about convergency of the infinite substitutions, i.e., the r-function is considered to be a 
formal series in A's (in the left hand sides, A's are implicitly present in the definition of discrete flows). 
In known examples of algebro-geometric solutions, the T-function is a true function, not merely a formal 
series. In this case there are some additional restrictions to the domains of all variables and labels. 
They ensure the convergency of the infinite substitutions. Meanwhile, for algebro-geometric solutions 
elements of the graph Q acquire a transparent interpretation on a Riemann surface as punctures and 
cuts. Furthermore, the discrete time variables discribe discontinuities of the Baker- Akhiezer function on 
the cuts. 



We refer to the discrete flows X a Xp as elementary discrete flows. One may introduce more complicated 
flows which can be thought of as "superpositions" of the elementary ones. Specifically, fix several 
elementary flows, say, h,fa, ■ ■ ■ ,1m (here h = for some cti,[3i) and consider the r-function as a 

function of a new variable y as follows: r[y] = t(Ii + y, I2 + y, • ■ ■ , Im + y)- In the time evolution with 
respect to the new variable y the "elementary" variables k simultaneously get shifted by y while the 
others are constants. Let us call flows of this type composite discrete flows. 

To put it differently, let di be the vector field corresponding to the flow U. Then the vector field 
corresponding to the composite flow y is d y := di, so 

exp(9 2/ )r({/ Qift }) = T({l aij3i + 1}) exp(<9 y ) . 



However, one should be careful since due to ( |3.2| ) the simultaneous shift of l a p and lp 1 is equivalent to 
an elementary flow. 

The precise definition is as follows: 
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• Composite discrete flows are labeled by finite sets of vectors {X ai X/3 i }, i = 1, 2, . . . , M such that 
0i 7^ ctj for any Let y be the corresponding time variable, then the evolution is defined by 

T M = r ({^/3 l + y}), 
where la.^ and other elementary variables are supposed to be constants. 

The distinction between elementary and composite flows can be extended to the continuous flows, 
too. For the reason which will be more clear later, it is natural to consider the continuous times t^ 
as elementary flows. At this stage we motivate this definition by the fact that due to Miwa's rule 



(3.3) they can be obtained as the result of a scaling limit from discrete elementary flows. Similarly, 
higher continuous times ty' with j >2 are limits of composite discrete flows. Therefore, we call them 
composite. 

Let us summarize. We introduced several notions and definitions which are extensively used through- 
out the paper. First of all, a partial classification of flows and time variables has been suggested. We 
have defined discrete and continuous flows and distinguished between elementary and composite flows. 
To any particular equation one may assigne a graph which explicitly shows the kinematical structure of 
the equation and possible constraints imposed on the flows. In order to make this clear, we give some 
examples. 

3.2 Examples 

Here we illustrate the above notions by familiar examples. To make the graphs of flows more informative, 
let us add a new graphical element: fat dots mean that the corresponding continuous times are non-zero^. 

In the KP hierarchy the graph of flows consists of one "fat" point with corersponding continuous 
"times" {tj}. The set of discrete flows is empty. The r-function is r(t) = r(ti, t2, ■ ■ ■)■ 

In the 2DTL hierarchy the graph of flows consists of two " fat" points with the corresponding times 
{tj} and {tj}. The discrete flow associated to the vector connecting the two points is the discrete "time" 
n of the 2DTL r-function T n (t; t) . 
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KP 2D Toda 

The graphs of flows for the discrete KP and 2DTL equations are as follows: 




Discrete KP (1-st eq.) Discrete 2D Toda (1-st eq.) 



All continuous times are swithed off. In both cases only three independent discrete flows are swithed on. 
This agrees with the continuous case where first non-trivial equations of the KP and 2DTL hierarchies 

4 Taking into account Miwa's rule (iii). it, w o uld he more precise to say that the continuous times can not be made equal 
to zero by a transformation of the form (|3.3|), (tj 
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have three independent variables ii,i2,^3 and t\ 1 t\, n, respectively. In the continuum limit, all the lines 
except the vertical one in the 2DTL figure shrink up to fat points. 

The discretized KP and 2DTL hierarchies are illustrated by the following graphs: 






Discrete KP (the hierarchy) 



Discrete 2D Toda (the hierarchy) 



Higher equations of the hierarchies involve more than three elementary flows. Labels of these flows are 
analogous to the number of a higher flow in the continuous hierarchies. The labels are complex numbers. 
This looks like a kind of Fourier duality between a parameter marking equations of the hierarchy and 
the time variable corresponding to a particular flow: continuous flows are marked by a discrete "label" 
whereas discrete flows are marked by a continuous label. 



4 Discrete Zakharov-Shabat representation of Hirota's equa- 
tion 



The reformulation of classical nonlinear integrable equations as flatness conditions for a two-dimensional 
connection is the basic constituent of the theory. The flatness means that subsequent shifts along any 
pair of the time flows commute. These conditions are known as Zakharov-Shabat equations or zero 
curvature representati on. In the paper |lj R.Hirota gave an example of the discretized zero curvature 



representation for eq. (IT). In the physical language, the discrete connection is a lattice gauge field. The 
approach emphasizing the relation to gauge field theories on the lattice was developed by S.Saito and 
N.Saitoh pi]]. We present these results in a modified form which makes the theory completely parallel 
to the 2DTL theory §§. 

The discrete zero curvature condition is equivalent to commutativity of certain multivariable differ- 
ence operators. The existence of such a "commutation representation" is a hall-mark of itegrability. At 
the same time, if a commutation representation exists, it is not unique. In particular, there are different 
(in fact infinitely many) ways to represent HBDE as a zero curvature condition. 

The general scheme is as follows. Choose any time flow to be the "reference" one, i.e., the one in 
which all the M-operators are going to act as differential or difference operators. Commutativity of 
the flows means that any pair of such M-operators obeys a compatibility condition which is just one 
of the Zakharov-Shabat equations. This fact allows one to relate different hierarchies to each other. 
In general M-operators are pseudo-difference^] or difference operators with matrix coefficients. Here we 
consider only the case of difference operators. More general examples are given later in Sect. 6. When the 
reference flow is taken to be an elementary one, the coefficients are scalar functions. Sect. 4.3 contains 
an example of the zero curvature condition for HBDE realized by 2 x 2-matrix difference operators. 

5 We use this shorter name for what is usually called "quantum pseudo-differential operator". 
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4.1 Basic M-operators 

So far all elementary discrete flows enjoyed equal rights. No one of them was better than any other 
one. We are going to break this equality of rights and distinguish a reference flow. It may be any flow 
including composite and continuous ones. For simplicity, we start with the case when the reference flow 
is discrete and elementary. Other cases are discussed later. 

The idea is to assign difference operators to flows. These operators act to functions of the reference 
flow variable. We call them M-operators. In this section we consider the simplest M-operators which 
are basic blocks of more general ones. 

Let us specify the notation and take the reference flow to be A0A1. The double index notation is 
inconvenient for practical purposes. Dealing with a limited number of Hows, it is worthwhile to give 
them simpler though less systematic names. Unless otherwise stated, the letter u will be reserved for 
the reference variable corresponding to an elementary discrete flow. So we set 

u = l i. 

Let A2 be any label different from Ao, Ai. In this situation, the graph of flows is a triangle: 




Its sides A0A2 and A1A2 define flows which we call adjacent to the reference flow u in the obvious sense. 
In general, a flow A«A 7 (resp., A^Ay) is said to be left adjacent (resp., right adjacent) to the flow AqA^. 

Coming back to the triangle graph of flows, we set Z02 = I, W2 = V • The r-function will be denoted 
by . There are only two independent flows. According to (ii), we have 

T%t X =r 1 ^ 1 '. (4.1) 

Coefficients of M-operators are expressed through r. Let us take u, I as independent variables. By 
definition, the M-operator assigned to the left adjacent flow I is 

T i T i+i 

M u = e A 2 l+l 1 < ( 4 - 2 ) 

T « T u+1 

where the coefficient A^ 1 is expressed through the three labels: 

A" 1 = 7— — j- - T~~~r' A ? 2 = - A ° 1 - ( 4 - 3 ) 

A2 — Ao Ai — Ao 

The shift operator e du has standard commutation relations with functions of u: e ±a "/(u) = f(u ± 



l)e ±9u . Note that M™ = e du . It is implied that the r-function in (4.2) might depend on all the other 
variables which are swithed off in this particular case. When they are swithed on, they enter eq. (4.2) 
as parameters. Their values are the same for each of the four r-functions in the ratio. As a rule, we do 
not indicate them explicitly if this does not lead to a confusion. 

Once the Af-operator for the left adjacent flow is written, it can be translated into the one for the 



right adjacent flow I' by passing to independent variables u, V . In eq. (4.2), the (implicit) argument V is 
the same in each r-function. Using ( |4.1| ), we rewrite them in such a way that I is the same and implicit. 
The rule (ii) tells us that the shift V — > Z' + l is equivalent the simultaneous shifts u — * u+l and / — > l+l. 
Thinking about M-operators as generating shifts of the discrete variables by 1, it is natural to define 
the M-operator for the right adjacent flow V as follows: 

Ml' = e~ 9 "Mi, 
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or, more explicitly, 

tV _ 1 \oi 'u+i 'u-i _a 



It is also useful to introduce 

M l u = e- d 'M l u , M l i=e- 9 <Mi (4.5) 

which are difference operators in two variables. It trivially follows from the construction that they 
commute: 

[MlMt]=0. (4.6) 

We have defined M-operators for elementary discrete flows adjacent to the reference one. In this 
case they have especially simple form. They are first order difference operators in u. M-operators 
corresponding to composite and non-adjacent flows have more complicated structure. 



Let us comment on continuous reference flows. According to (3.3), the continuum limit in u means 
e a " — > 1 — A _1 9 fl + 0(A~ 2 ), where t\ is the first continuous flow with labels Ao and A = (Ai — A ) _1 . 
The limiting form of the M-operator (4.2) as A — + oo is 



i+i 



M*'> = d tl - d tl log — -j (A 2 - A 



i-i 



This is a first order differential operator in the reference continuous time veriable t\. It generates shifts 
in the discrete variable I. The fact that it is an operator of the first order suggests to call the continuous 
flow ti elementary (see the end of Sect. 3.1). 



4.2 Discrete Zakharov-Shabat equations 

Two independent flows are not enough for deriving bilinear equations. Non-trivial bilinear equations for 
r arise starting from 3 independent discrete flows, in which case the graph of flows should contain at 
least 4 vertices. So let us fix 4 labels A Q , a — 0, 1, 2, 3 and consider the general graph with 4 vertices: 




The simplified notation for the flows is clear from the picture. Like in (4.3), we set 

1 1 



A «/3 



A^ — A Q ' 

for all possible values of the pairwise distinct indices. 



A /3 



7 



(4.7) 



Let AqAi be the reference flow, as before. The left (resp., right) adjacent flows are A0A2, A0A3 (resp., 



A1A2, A1A3). Each of them has its own M-operator of the form (4.2) (resp., (4.4)). 
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The key point is to extend the trivial commutation (4.6) to all the flows in the graph adjacent to u: 

[Ml AC] = [Ml M{] = [Ml M™] = . (4.8) 



In contrast to eq. (4.6), these are non-trivial requirements which give bilinear equations for r. Written 
in terms of M-operators, commutation relations (4.8) are discrete Zakharov-Shabat equations. 

In the following proposition we use the notation like M l u = M l u (u, 1,1,.. .) to indicate the dependence 
of M-operators on the variables. 



Proposition 4.1 The discrete Zakharov-Shabat equations 

M™(m, I + 1)M< (m, I) = M l u {m + 1, l)M™(m, I) . 
M™(m,T+ l)M T u (m, I) = M l u (m + 1, l)M™(fh, I) . 

mIu + i, i) = l ) M % 

are equivalent to the following bilinear relations for t: 



\01_2+l,m i,m+l 
A 2 T u T u+1 



A 3 V T u+1 



\01 „J,m+l l+l,m 
A 3 ~u T u+1 



\01_J,m+l l+l,m 
A 2 T u ~u+l 



Hi(l,m;u)T u 



l+l,m+l l,m 
T u+1 



l,m Z + l,m+l 



■H 2 {l,m;u)T^T? + \ 



0. 



0, 



,01 1,1 + 1 1+1,1 
K 2 T u+1 T u-1 

respectively, where Hi are arbitrary functions such that Hi(l,m;u + 1) = Hi(l,m;u). 



\0i 1,1+1 1+1,1 \0l J,l+l J+l,l , tt a J. „\ J,(J+l,f+i n 

A 3 T u r u -*2 T u+l T u-l + H 3{IJ;U)T^ T u =0. 



(4.9) 

(4.10) 
(4.11) 

(4.12) 

(4.13) 
(4.14) 



Proof. The proof consists in the straightforward commutation of M-operators. The M-operators read: 



Ml = e°" 



T i J+i 

v 01 'u'u+l 
^3 l+l "J ' 

Ttl T u+1 



m,'< 



m _m+l 

,oi 'u-i'ti+1 „-g u 

A 3 m+l ^ e 



(4.15) 



(4.16) 



M™ is given by ( 4.15 ) with the cha nges I — > m and 3 — ► 2, M Jt is given by ( 4.16 ) with the changes m — > Z 



and 3^2. The details for eq. (4.S) are given below. 
Eq. (fh|) reads: 



-m,l+im+ld + l 

\01 U M+l 

A 2 m+1,2 + 1 m,2+l 
T « T u+1 



_m+l,im+l.'+l 
\01 « 



'u+1 



m+l,i+l m+l, I 



'u+1 




r m,lm,l+l ' 
'u 'u+1 
m.l+1 m,l 
r « ' T u+1 , 

r m,lm + l,l ' 
'u 'u+1 
2 m+l, 2 m,l 
Tu T u +1 , 



(4.17) 



The terms e 2du and those which do not contain shift operators cancel automatically. The comparison of 
the coefficients in front of e du yields 



m,l+l 



-i,i ' 



T u + 1 I x 01 T u+2 _ 01 T u+2 
m,l I A 3 m ,l+l 2 m+l,l 



'u+2 \ 

m+l,i + l 
T u+1 
m+l,l+l 



'u+1 



'u+1 



i01 'u 



r m+l,Z 



\01 T u 



m,l + l 



m+l,l 



m,l + l 



'u+1 



'u+1 



or 



\01_m+l,Z m,Z+l 
A 3 T u T u+1 



m,l m +l,l+l 

T u+l T u 

\01 m+l,/ m,L+l \ oi m,/+l m+l,/ m,Z m+l,Z+l 

A ^ T"„,_L1 T., i n An T,, i 1 i n T_, i oT,, i 1 



01_m,;+l m+l,i 
'u+1 



A 2* T < 



(4.18) 



(4.19) 



*3 'u+1 'u+2 ^2 'u+1 'u+2 'u+2 'u+1 

The denominators in both sides differ from the numerators by the shift u — *■ u+- 1. Therefore, their ratio 
is a " quasiconstant" in u, so the equation is equivalent to eq. (4.12). This completes the proof. 
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Now we restore the "equality of rights" of elementary flows by imposing the requirement that 
Zakharov-Shabat equations should hold for any choice of the reference flow. For example, let I be 
the reference flow. Construct .M-operators for the flows u,u adjacent to I (see the figure). Then we 
require that the operators Mf, A4" 1 Mf, M~[ m commute with each other (of course some of them 
commute automatically due to ( |4.6| )). Note, however, that .M-operators constructed with respect to 
different reference flows are not required to be commuting, e.g. [M™, M™\ ^ 0. 



Theorem 4.1 Let x be any one of the elementary flows shown in the picture and let v, v be the corre- 
sponding left and right adjacent flows such that x,v,v are independent. Then the commutativity condi- 
tions 

[M%, M%] = 

imposed simultaneously for any three independent reference flows x are equivalent to the equations 

\03_J,fn+l_J+l,m , \ 01l+l.ml,m+l , \ 02J+l.m+l l,m _ n (A IDA 

A l T u T u+1 + A 2 T u ' T u+1 + A 3 T u ' T u+1 - U ) 

A l T u T u+1 + A 2 T u T u+1 + A 3 T u ' T u+1 = > I 4 - 21 ) 

\ 01^1,1+1^.1+1 J \02i,Ii+i,i+i a,oi i,T+i i+iJ (a on\ 

A 3 T u T u - A 3 T u T u = A 2 T u+1 T u-1 ■ ( 4 - 22 ) 



Sketch of proof. By virtue of the previous proposition, it is enough to show that the functions Hi are 
constants: Hi = -H 2 = -H 3 = Xf. This can be done straightforwardly by writing down bilinear 
equations arising from Zakharov-Shabat equations for M-operators corresponding to each choice of the 
reference flow and demanding their consistency with each other. 

We see that eq uatio ns ( [4.20[ ), fl4.22|) coincide with the KP and Todadike forms of HBDE ( |2.1l| ), 
respectively. (Eq. ( |4.21 ) coincides with the KP-like form after the change u — ► —u.) T he thr e e equ ations 
differ by the choice of the independent variables only that agrees with substitutions ( [2.14 )-( p.l9| ). The 
transition from one triple of independent variables to another should be done according to the rules (i) 
and (ii) (see (3.1), (^^)). Using these rules, it is easy to see that the three equations (4.20)-( I.2l| ) arc 
equivalent to each other. 

The 4-variable MKP-like form ( [2.10 ) of the Hirota equation follows from (4.20) by applying the 
rule ( |3.2| ). Namely, fix an extra label and consider the flows /zoA Q , a = 0, ... 3 with time variables 
p a . From (|3~2 ) we have, for instance, t™^ 1 ^ 2 = t™^ 2+1 and so on. This change of variables converts 
equation ( |4.2C ) into eq. ( |2.10[ ). 



4.3 Matrix realization of the zero curvature condition 

We restrict ourselves by giving an example which illustrate the general scheme outlined in the introduc- 
tion to this section. 

Consider the graph of flows which is a reduced version of the one from Sect. 4.2: 

Ai ; A 2 



Ao 
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The variables which we do not need here are swithed off. The simplified ad hoc notation is clear 
from the figure. This choice of independent variables corresponds to the discretized 2DTL equation. 

Our goal here is to write the zero curvature condition with another choice of the reference flow. 

> > 

Specifically, let it be the composite flow labeled by the pair of vectors A0A3, A1A2. Let y be the cor- 
responding "composite" time variable. According to the definition given in Sect. 3.1, the r-function 
depends on y as follows: 



In other words, we set, by definition, d y := di + di, so the shift operator e dy acts to the r-function by 
shifting 1,1 simultaneously: 



e dy T i,i = T l+l,l+x e & y 



Introduce the following difference operators with 2 x 2-matrix coefficients: 



/ 



L n (i,T) 



+ v 



<n 'n+1 
1+1,1 1,1 
r « T n+1 



Z+l.Z + l 



/./,- 



'n+1 



'n— 1 



.1+1,1+1 1+1,1 



-[IV- 



Z+l.Z+l 
r n+l 



J.I 



1,1 
T n+1 



M n {l,l) 



1,1+1 

r n+l 
1,1+1 



V 



Z+l.Z+l 
r n-l 
1+1,1 + 1 



1,1+1 Z+l.Z+l 
T n+1 T n-1 

1 1,1 + 1 1+1,1+1 / 
Tn 7~n ' 



where we put /i = A!] 1 , v = A3 1 for the sake of brevity. 

Proposition 4.2 The matrix discrete Zakharov-Shabat equation 

L n (l,l + l)M n (Z,0 = M n+1 (l,l)L n (l,l) 
is equivalent to the bilinear relation 



.Z+l.Z+l 



— H n (l, T n r i 



1,1-1+1,1+1 _ 



j,i+i 1+1,1 

n+1 T n-1 1 



where H n (l,l) is periodic in n with period 1: H n+ i(l,l) = H n (l,l). 



(4.23) 



(4.24) 



(4.25) 



This bilinear equation coincides with eq. ( 4.14 ). We omit the proof since it is absolutely straightforward 
after the L-M pair is given. A way to derive matrix M-operators from the scalar ones is discussed in 
Sect. 5. 

Like in Theorem 4.1, the validity of the zero curvature condition for M-operators constructed with 
respect to all possible independent reference flows implies the bilinear equation with a fixed constant 
function H . It has the form ( 2.11 ). 

Remark 4.1 In the 2DTL interpretation, the operator M n generates evolution in the chiral discrete 
"space-time" variable / whereas L n generates shifts along the n-lattice. In our scheme, both M n and L n 
are " Af -operators" rather than " L-operators" . We write L n according to the tradition which is justified 
by the case when an additional reduction of the 2DTL is implied. 

It is instructive to look at the continuous version of this zero curvature condition. It provides the 
zero curvature representation of the 2DTL with the composite continuous reference flow defined by the 
vector field d y := d tl (see Sect. 3). This representation naturally arises when one embeds the 2DTL 
into the 2-component KP hierarchy. The Zakharov-Shabat equation 



dtiL n — M n +xL n — L n M n 



(4.26) 



14 



where 



T n +1 



M„ 



T n -1 



r n +l 



(4.27) 



is equivalent to 



d^dt, log- 



'ra+1 



(4.28) 



which is the 2DTL equation in the bilinear form. 

Remark 4.2 The ID Toda chain (1DTC) is a reduction of the 2DTL such that t„ does not depend 
on ti + ti, i.e., d y commutes with r n . Therefore, in this cas e d v can be considered as a c- number. 
Identifying it with the spectral parameter, one recognizes eqs. ( 4.27 ) as the standard L-Af-pair for the 
1DTC realized by 2 x 2-matrices depending on the spectral parameter (see e.g. f^|). 



5 Linearization of the Hirota equations 

Zero curvature conditions studied in the previous section are equivalent to compatibility of an overde- 
termined system of linear difference equations for a "wave function" ip. These linear equations are called 
auxiliary linear problems (ALP). They play a very important role in the theory. Common solutions to 
ALP carry the complete information about solutions to the nonlinear equations. All the properties of 
the latter can be translated into the language of the ALP. This is what we mean by the linearization of 
HBDE. 

In accoradance with the diversity of zero curvature representations there are many types of the ALP. 
This section deals with the most important examples. 



We begin with the scalar linear problems associated with the M-operators (4.15), (4.16) for el- 
ementary discrete flows adjacent to the reference one. They are simple first order linear difference 
equations with coefficients expressed through the r-function. The formal solution of a special form is 
called (formal) Baker- Akhiezer function. It depends on a spectral parameter. Baker- Akhiezer functions 
are formal analogues of Bloch solutions. The formula for Baker-Akhiezer functions in terms of the r- 
function was suggested for the first time in |^4|. General solutions to the ALP are linear combinations 
of Baker-Akhiezer functions with different spectral parameters. In a similar way, one may define dual 
Baker-Akhiezer functions as formal solutions to the linear problems for adjoint operators. 

Given a solution to the ALP, one may consider Backlund transformations. Furthermore, a "duality" 
between coefficient functions and solutions of the ALP allows one to define a chain of successive Backlund 
transformations described by the Backlund flow. We consider two types of the Backlund flows. It is 
shown that in the particular case when the solutions to the ALP are Baker-Akhiezer functions the 
Backlund flows can be identified with elementary discrete flows adjacent to the reference flow. 

There is a "gauge freedom" in the ALP which can be fixed by certain normalization of ip. Mostly 
we use the gauge which leads to the simplest possible form of the linear equations. Another choice - 
zo-gauge - is briefly discussed in Sect. 5.4. This gauge makes equations more symmetric for the price of 
introducing an auxiliary point zq G C and complication of coefficient functions. 

The ALP associated with the matrix M-operators are also discussed. In fact matrix M-operators 
can be most conveniently derived using the scalar ALP. The matrix linear problems are obtained by 
rearranging the scalar ones. More precisely, in order to rearrange the scalar ALP in such a way that the 
reference flow is taken to be composite, one has to pass to difference operators with matrix coefficients. 
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5.1 Scalar linear problems 



The commutativity of .A/f-operators (4.8) implies that they have a common set of eigenfunctions. Equiv- 
alently, the discrete Zakharov-Shabat equations (4.9)-( 4.11 ) for M-operators imply compatibility of the 
linear problems 

M 1 u tP 1 < T (u)=iP 1+x >\u), (5.1) 

MV''(u)=?i a+1 (ti) (5.2) 

for any elementary discrete flows I, I adjacent to u. Note that the "eigen value s" put equal to 1 in the r.h.s. 
can be made arbitrary by changing normalization of ip. Our choice in (5.1), (5.2) is most convenient in 
the purely discrete case though it does not lead to a smooth continuum limit. 

More explicitly, eqs. flsTil ), Q read (see (|1|), (^16|)): 



^{u + 1) - A° 1 V l ' T (u)i> lJ {u) = ^ +U (u) , 
VM(u) - \fC l \u)i, l \u - 1) = iP lJ+1 (u) , 



where 



V M (u) 



C u (u) 



-J,L- i+1 >' 

T u T u+1 

1+1,1 1,1 

r « T u+1 

1,1+1 1,1 

T u+1 T u-1 

1,1+1 1,1 

Tu T~ u 



These formulas become more symmetric in terms of the " unnormalized" wave function 



J,i 



i.i 



Plugging this into (|5.3| ), ( p.4| ), we get: 



l,l J+i,i 



l+i,l 1,1 iOi 1+1,1 1,1 _ ■ . 

T u Pu+1 A 3 T u+1 Pu — T u+lP- 



M+iJ,i \0i 1,1+1 1,1 _M n l,l+l 

T u Pu ~ A 2 T u+1 Pu-1 — T u Pu 



(5.3) 
(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 
(5.9) 



Let us show that a carefully performed continuum limit of these equations yields the familiar linear 
problems for the 2DTL. Let 

A3 = A + e, A 2 = Ai + e, e, e -> , 



so 



1 - ec» tl , 



1 



edf. 



by virtue of Miwa's rule (3.3) applied to the discrete flows A0A3, A1A2, respectively. Let us change the 
normalization of the ^-function introducing the (^-function as follows: 



^(u) = (A - A 3 )V <<(«) 

so the linear problems read 

(A - X 3 )^ T (u + 1) - A(A 3 - \i)V l > T (u)<p l > T (u) = v l+1 '\u) , 

(Ao - X 2 )<p l >\u) - A(A 2 - AiJC'fyyfy - 1) = (A - X 2 )^ I+1 (u) . 

where A = (Ai — A ) _1 . We get as e, e — > 0: 

T u+1 1 , , 

ip{u) 



dtMu) =ip(u+l)+[ X + d tl lo. 



(5.10) 



(5-11) 



(5.12) 
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The transformation 



X -u e -Xuti 



eliminates the constant A turning the linear problems into the familiar form 

d tl ip{u) = ip(u + 1) + v(u)tp(u) 

d^fiu) = c(u)(p(u - 1) . 

Here 



(5.13) 



v(u) = d tl \o\ 



T u +1 



c(u) = 



T u +lT u -l 



(for the sake of notational simplicity we use here the same letter for the transformed function). The 
Zakharov-Shabat equation 



[d tl - e a * 

yields the 2DTL equation in the form 

d tx df x log 



v{u), d h -c{u)e-^\ =0 



T u +1 T u+ \T u ^i T U+2 T U 



'u+1 



(5.14) 



These are continuous analogues of eqs. fl4.1l|) , ( |4.19| ), respectively. 



5.2 Backlund transformations 



The ALP in the form (5^8), (5^) have a remarkable property emphasized in ref. |2l|| : they are symmetric 
under interchanging r and p. Furthermore, one may treat them as linear problems for the function r, 
the compatibility condition being a bilinear equation for p. This equation is again HBDE of the same 
form. In ref. (2lj, this fact was refered to as the duality between "potentials" r and "wave functions" p. 
This "duality" shows up most transparently in the fully discrctized case. 



More precisely, rewriting eqs. (5J3), (5A) as linear equations for 



^ M H=^ 1 ^=(^ +ij+i («+ir J 

Pu+i 



(see (p.7|)), we get 



\fV l \u))^\u)=^-^{u), 



71.1 



1-1,1, 



1 - \fC l '\u + l)e°- ) Tp l ' l (u) = ^' l - l {u) 



7.1,1-1 1 



(5.15) 
(5.16) 



where V and C are given by the same fo rmula s (5.5), ( p.q ) with p in place of r. The difference operators 
in the l.h.s. are adjoint to the operators ( 4.15 ), ( 4.16| ) with t — > p. The formal adjoint operator is defined 
by the rule (f(u)e kdu y = e~ k9u f(u). It then follows that the compatibility conditions are described by 
Theorem 4.1 with r replaced by p. 

Therefore, passing from a given solution r to p we have got a new solution to HBDE. This is a 
Backlund- type transformation also known under the names "Darboux" or " Backlund-Darboux" trans- 



formation. For a comprehensive discussion of transformations of this kind see |25|. The bilinear form of 
Backlund transformations was suggested by R.Hirota ||. 

One may repeat the procedure once again starting from p and, moreover, consider a chain of successive 
transformations of this type. It is natural to introduce an additional discrete variable b to mark steps 
of the "flow" along this chain and let r^ l b , p l ^ b be r and p at 6-th step, respectively. The first Backlund 
flow is defined by 

T u!b+i=Pu,b- ( 5 - 17 ) 

This mean s th at rat the next step of the " Backlund time" b is put equal to a solution p of the linear 
equations (5^8), (5J5). Then these linear problems become bilinear equations for t&: 



T u,b T tt+l,b+l 



\01 l+l J _ 1 l+l 

A 3 T u+l,b T u,b+l T u+l,b T u ,b+l 



0. 



(5.18) 
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J+lJ 



b T u,b+i( u ) - t uAu)t u \ +1 = A 



\01i + l I 



-l,b' u-l,i>+l ' 



(5.19) 



where I (resp., I) in eq. (5.18) (resp., ( 5.19| )) is skipped. 

Similarly, defining the second Backlund flow (the Backlund "time" is now denoted by b), 



1,1 

T u,b+1 



l,l 



we get from (5.8), (|5.9|): 



\01 J+i / _ 
A 3 T u h T u,b+1 



u,b' u ,b+l 'u+1,6+1 'u-1,6 : 



u,b u+l,6+l 



J + 1J 



-1,6+1 + A 2 T „+i S T «,6+1 - U 



(5.20) 

(5.21) 
(5.22) 



In these equations one immediately recognizes different forms of HBDE. A time discretization of the 
Toda chain by means of Darboux transformations was considered in p(i|| . 

The Backlund flows can be defined by a zero curvature condition. Given any solution ij> to the ALP 



(5.3), (5.4), introduce the operator 



B J 



ip(u) 



(5.23) 



Then eq. (5.18) is represented as the commutativity condition A4 U ] = 0. A similar S-operator exists 
for the second Backlund flow. 



5.3 Baker- Akhiezer functions 



Each of the ALP (^j), ( p.4[ ) is a first order linear difference equation in two variables. Assuming HBDE 
(4.20)-(4.22) hold, we are going to construct a 1-parametric family of their common solutions of a special 
form. These solutions ip{ u ) = i ) { v <\ z ) are called Baker- Akhiezer functions. They depend on the spectral 
parameter z € C. 

Let us swith on an extra elementary flow shown by the dotted line: 



Ai 



A 



The corresponding time variable is p z . Let 



A/? — A 



Note the identity 



A?' 3 A? 7 = A£ 7 A? 7 . 



(5.24) 



(5.25) 



Then, assuming the 3-term Hirota equations hold for for the triples (it, l,p z ) and (it, l,p z ) of independent 
variables, 



^(u;*) = (A-rm 



\02 \ 1 T- 1,1 
01\uf\03\l I A z \ 'u,p z + l 



A" 1 



i,i 



(5.26) 



Pz=0 



is a formal co mmon solution to eqs. (5.3), (5.4) for any z. Indeed, under th e sub stitution ( 5.26 ) eq. (5.3) 
turns into eq. ( 4.20 ) for the triple (u,l,p z ) while eq. (5-4) turns into eq. (122) for the triple (u,l,p z ). 
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Therefore, the new label z is identified with the spectral parameter. Formula (5.26) for the ^-function 
coincides with the Japanese formula jl5), (24) because due to ( |3.3| ) we have 

1,1 

^~u,p z -\-l 



U 



r£''(0) 



Pz=0 



The general solution to the ALP can be represented in the form 

ip(u) = / d 2 zfj,(z)ijj(u; z) 



(5.27) 



with arbitrary "measure" /i(z) on the complex plane. In other words, this is a linear combination of 
Baker- Akhiezer functions with different spectral parameters. 



Note that the " £>-operator" ( 5.23 ) in which ip is taken to be the Baker-Akhie zcr fu nction coincides 
with an .M-operator. Indeed, we have the following formula for the M-operator ( 4.15 ) in terms of the 
Baker-Akhiezer function: 

Ml = lim - 



(5.28) 



The dual Baker-Akhiezer function ip* is defined by 



-i i,i 



U,p z —1 



1,1 



(5.29) 



Pz=0 



It satisfies the equations 

(Ml(u -1,1- l)) f r\u; z) = z) , 

(m T u (u - 1, 1- lj) tp* T (u; z) = ip^iu; z) , (5.30) 
where the difference operators in the l.h.s. are formally adjoint to the operators ( 4.15| ), ( |4.16[ ). 



5.4 The 2 -gauge 



Equations ( p.3[ ), ( p.4| ) imply a specific choice of normalization of the -i/>-function. Indeed, multiplying tp 
by any function, one can change the form of the equations. This is a kind of a "gauge freedom". There 
is no canonical way to fix the gauge. The gauge that we systematically use throughout the paper is the 
most economic one in the sense that the ALP for discrete flows have the simplest possible form. Here 
we are going to discuss another choice which has its own advantages. 

This more general gauge requires to fix an additional point zq G C different from the vertices of the 
graph of flows. The gauge is defined by the following normalization condition for the Baker-Akhiezer 
function ^(u; z): 

V{u;z ) = l. (5.31) 
We call it zq- gauge. Given this condition, it is natural to represent \E' in the form 

*(«;»)- VJ(U;Z) 



i>(u; Zq) 

and rewrite the ALP ( |5.3| ), (|]J) for ip in terms of ^. In this way, we get 

(A°Jtf(u, 0e S " - X° z lW(u, l)e 9 ') *(u; z) = \f*(u; z) , (5.32) 

where 

T i+i T i r '+i r l 

U( u> l) = ».P0 T «+1,P0 + 1 W ( /} = u,po + l «+l,Po _ (5 _ 33) 

T T T T 

'«,PO+l U+l,po 'u,po + l M + l.PO 

A general prescription for writing down the equation should be clear from comparison with the picture. 
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With this prescription at hand, a similar equation can be written for any pair of flows such that one of 
them is reference and another one is left adjacent. 

A nice feature of the z -gauge is that there is no need to care about equations for right adjacent 
flows. They are automatically produced by the same prescription if one changes the "orientation" of 
the reference flow (i.e., consider A1A0 as the reference flow). In order to express everything in the same 
variables, one should appl y th e rules (3.1) (it — » —u) and ( [3^2] ) (po — * Po)- We stress that eq. (5^4) can 
not be obtained from eq. ( |5.3|) in this way. In that gauge we need two types of equations separately for 
left and right adjacent flows. 

The Baker- Akhiezer function that solves all these equations in the Zo-gauge has the form 



*(«;*) = 11 

a/3 



A «/3 



'P0,Pa; + l 



(5.34) 



PO =Pz=0 



As usual, the r-functions depen d on a ll th e skipped variables as parameters. Due to (5.25) the form of 
the prefactor is consistent with (3.1), (3.2). 

Our previous gauge is a limiting case of the zo-gauge as zq — > Ao- However, the limit is singular and, 
therefore, it needs a rcgularization. As a result, the symmetry of the zp-gauge becomes broken. 



5.5 Two-component formalism 



Linear equations ( p.3[ ), (5.4) can be brought into another form in which they become first order partial 
difference equations for a 2-component vector function. Their compatibility yields matrix Zakharov- 
Shabat equations presented in Sect. 4.3. 

Let us use the notation of Sect. 4.3 and, accordingly, denote i\) n = ip(n), V^' 1 = V Ll (n), = C l ' l (n). 
Then eqs. (Q, ([[J) read: 



These equations allow us to find out how the vector 

ill 

transforms under shifts in n and I. 

Combining eqs. ( 5.35] ), ( 5.36 ) we have, for instance (here d y = 



di+di): 



i i,i 



= (e ^ + vVk 1 + ~ »vC l +^V^ l li 



(5.35) 
(5.36) 

(5.37) 



(5.38) 
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Proceeding in the same way, we get: 

1 

1 -fxCk 1 



i i,i 

Vn+1 

^ij+i 



,1,1+1 




(5.39) 



The operators in the r.h.s. provide a matrix L-M-pair for HBDE which differs from (4.23) by a "gauge" 
transformation. Recall that the Baker- Akhiezer function ( [5.26| ) has r n in the denominator, so the two 
components of the vector (5.37) have different denominators. In the 2-component formalism, it is natural 
to require the denominators to be the same. This condition partially fixes the gauge. 



Thus, introducing the vector (ip n , Xn) with the second component 

1,1 

n — l 1 1,1 



T, 



we rewrite eqs. (5.39) in the form 



L n (l,T) 



M n (l,l) 



^n+l 

1,1 
Xn+1 



^1,1+1 
1,1+1 

An 



(5.40) 

(5.41) 
(5.42) 



with the L and M-operators given by eqs. (|4.23j) . 

Eq uations ( 5.41 ), ( 5.42 ) imply some useful difference equations for . Excluding Xn, we get from 



^l+U+i 



i>n+l 



T n T n+1 T n-1 1 n+1 

V , , ; f , f h (Jt- 



1pli l + fiU 



1+1,1 1,1 
T ™ T n+1 



. ,-1,1 + 1 l+l,i 
Tn. Tn. 



1,1 1+1,1 + 1 
r n-l T n+l 



1,1 



1,1 1+1,1+1 ' 

Tn Tn 



(5.43) 



Similarly, excluding Xn from eq. (5.42) and using HBDE ( |2.1l| ), we get a 4-term equation for ip n 



4>- 



1+1,1+1 



1>- 



i+i,i 



r i,i+i 1+1,1+1 



M 1+1,1+1 

'n 1 n+1 



1+1,1 + 1 1,1+ 
r « T n+1 



I l.l + L i / \ n ■ n+l ; 



1.1,1 



(5.44) 



'n+1 



Here it is im plied that r„> satisfies HBDE ( 4.22 ) in the corresponding notation. For more information 
on eq. ( 5.44 ) see the next section. 



The continuous analogue of eq. (5.43) is 



(d tl + dtjfn = <^„+i + d tl log 



T n +1 



fn 



Tn-lTn+1 



-fn-1 



(5.45) 



This is obvious from (5.12) (we write ip n instead of ip(u) according to the present notation). It is a 
discrete non-stationary Schrodingcr equation. Eq. (5.44) in the continuum limit turns into 



StiAxVn - <k x 10; 



T n +1 



dt^n 



Tn-lTn+1 



(5.46) 



It is a continuous two-dimensional Schrodinger equation in magnetic field. Its quasiperiodic solutions 
were studied in the paper |27j by means of the algebro-geometric approach. The corresponding theory 
for discrete two-dimensional equations similar to Q5.44 ) was proposed in p8[ . 



21 



6 Pseudo-difference M-operators 



In this section we study the general form of Af-operators which satisfy conditions of the zero curvature 
form with M-operators for elementary discrete flows adjacent to the reference flow. 

Starting from the scalar ALP for a pair of left and right adjacent flows it is not difficult to find the 
M-operators for non- adjacent flows. Indeed, it is possible to exclude the reference flow from the pair of 
linear equations. Then the right adjacent flow turns out to be written in terms of the left adjacent one. 
Considering the latter as a new reference flow, one obtains a general M-operator for any elementary 
discrete flow in terms of any (elementary discrete) reference flow. In general these are pseudo-difference 
operators, i.e., they contain negative powers of first order difference operators. 

This construction can be extended to more general operators which generate some new flows commut- 
ing with the elementary discrete flows. We call them adjoint flows. The corresponding pseudo-difference 
operators are constructed in Sect. 6.2 with the help of two arbitrary independent solutions ip, V* to the 
ALP and the adjoint ALP. 



6.1 M-operators for arbitrary elementary discrete flows and the correspond- 
ing linear problems 



Rearranging eqs. ( p.3| ), ( p.4| ), it is possible to find out M-operators for arbitrary elementary discrete 
flows, not only for adjacent to the reference flow. The idea is to exclude shifts in u and then consider I 
as a new reference flow. 



From (5.3) we have 



V> M (u + 1) = 4> l+l - l {u) + X u 3 L V l ' l (u)^ l (u) 



- ,/,i+i,Z 



OltrU r.A.iJJ / 



Plugging this into (5.4) written in the form 

^' T (u + 1) - ^ I+1 (u + 1) = XfC l '\u + l)^'\u) , 
in place of ip l,l (u + 1) and ip l ' l+1 (u + 1), we get 

^ l+l - l+l {u) + X^V^iu^+^u) = ^ l+l - l {u) + (\fV lJ {u) - XfC lJ {u + l))^\u) . 



Using ( [4.221 ), we find 



A^V M (u)-A^C'>+l) 
XfV l ' l+l {u) - \fC l+l \u) = A 



01 'u 'u+1 

1+1,1 1,1 + 1 
r « T u+1 



(6.1) 



Note that the first equality follows already from eq. QL11 ) which is a weaker condition than (4.22). 
Therefore, ip obeys the following 4-term linear equation: 



if 



l+i,i+ii 



_z,r+i 1+1,1+1 

l, l + 1 > l (r,A — \03 'u 'u+1 



(u) - ^(u) = A f— \ T+1 ^ l+1 (u) + A^^__^>) 



: / 1+1,1+1 



i02 '« 'u+1 



(6.2) 



'u+1 



'u+1 



in which we recognize eq. (5.43) from Sect. 5.5. 

The relation ( |6.l|) allows us to rewrite it in the form 



where 



'u+l 



A/ ee e Sl + A 01 



3 ' 



Aj ee e d ~' - 1 



(6.3) 



(6.4) 
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This equation looks like a discrete 2-dimensional Laplace equation in a curved space. It can be formally 
rewritten as 

/ 1,1+1 



A i -+A° 1 ^±i T A r 1 ^l T ]^(u)=0, (6.5) 



or, finally, 



1,1+1 1 i+i.j 

>u 'u 



* ,J+1 <°»= y-^^rwT^j )*»■ < e - 6 > 

To avoid a confusion, we stress that the operator inside the brackets acts to the variable I whereas u 
enters as a parameter. This operator should be identified with the M-operator generating the flow I 
realized as a pseudo-difference operator in I. 

Ai I A 2 



Ao 



In other words, letting I to be the reference flow, we, therefore, have got an M-operator for the flow I 
which is not adjacent to I. 

In the limit when the points Ao, Ai merge, the flow I becomes left adjacent to I. Let us demonstrate 
how the corresponding difference M-operator is reproduced from the r.h.s. of eq. ( |6.6] ). Let Ai — Ao = e, 
e — > 0. We have 

X 01 

-^^^l + e^+AH+OCe 2 ), 



T u±l 



LI 



I + O(e) : 



1,1+1 1+1,1 

T u+1 T u-1 



1+1,1 1,1+1 

Tu Tu 



1 - 



I - 



T l,l T l+l,l + l 
r l+l,I~l,I+l 



+ 0{e 2 ) 



as e — ► (in the last line eq. (4.22) is used), so the naive limit of the r.h.s. ( |3.6] ) is zero. However, we 
should take into account the change of normalization which is implied when the former flow I turns into 
a left adjacent flow to /. This is achieved by ijr — > (— e) l i/> , so in the limit the right M-operator 



A 



li:S 



T l,l T l+l,l + l 

T l+i,T T l,T+i 



is reproduced. 

For illustrative purposes, let us give continuous analogues of the abov e for mulas. Rather than to 
perform the limit directly, it is much easier to use the continuous version (5.13) of the linear problems 
from the very beginning. Making use of eq. ( f5.14| ), we find the analogue of eq. ( 3.2): 



d^di^iu) - v{u)di x y{u) - c(u)(f(u) = . 
The analogues of eqs. (|6.3D, (|6 ,6[) read 



(6.7) 



T u -1 



-d^—^-d^ifiiu) = ip(u) 



Tu+l 



Or, 



T u +1 „-l1"ti-l 



(6.8) 
(6.9) 



respectively. 
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6.2 Adjoint flows 



Finally, we extend the above scheme to incorporate more general flows which we call adjoint. Let 

A? = 1 + wA-V , A, = e 9 ' - 1 (6.10) 

be a pseudo-difference operator, where w, w* are as yet arbitrary functions of all the time variables. 
By a we denote the time variable corresponding to the adjoint flow that we are going to define. In this 
section, the reference flow is I. The M-operator for an elementary discrete flow p (see the figure) is 

Mf = e 9 > - X P 4^I . K - Af . (6.11) 
T i+i T i 




Proposition 6.1 The commutativity condition 

[e- d "A?, e- d 'Mf] = (6.12) 
holds only if w and w* satisfy the linear equations 

(e*+\ p Vr +1 W=uwf +1 

(6.13) 



(e' 9 ' + X p V p ' a ) W ; p+1 ^ u;w* p 



i > 



whe 



and uj is an arbitrary constant. 



p,a p+l,a 

vr 



r p,a _ T l T l + l 



>.a p+l.a 



The proof is by straightforward computation. Equations ( 6.1 3| ) are necessary conditio ns for vanishing of 
the pseudo-difference part of the commutator. Here are main steps of the proof. Eq. ( |6.12| ), 

(> - x p vr +i ) (i + Arvo = (i + < +1 Ar^r 1 ) (e di - W) . 

can be rewritten as 

A P\ V l V l J+ W l W l-1 ~ W l + l W l 

- I ' m P \ V P ' a+1 m P , nn P+1 \ A^ 1 m* P 
~ [Wl+l ~ A P V l W l J A Z W l 

- wf +1 A- 1 (w^ 1 - X p V p ' a w* p+l - tow^ . (6.14) 

Since the l.h.s. does not contain negative powers of A/, the r.h.s. should be zero. This condition implies 
eqs. (3.13). 

The ALP for the left adjacent p to the reference flow I and its adjoint read (cf. (|5.3|) ) 

i>P(l + 1) - \ P Vf^ p {l) = 4> p+1 (l) 

(6.15) 

ip* p (l - 1) - ApV^V^O = ^* p ~ 1 (i) 
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Whence we identify 



1), 



where "0, "0* are arbitrary solutions to the linear problems ( |6.15[ ). 
The operator ( 6.10| ) acquires the form 

A? = l + ^ a+1 (l)AY l ip* a (l + 1) • 



(6.16) 



The commutativity condition (6.12) is equivalent to a nonlinear equation for rf ' a . The adjoint flow a 
is defined by two arbitrary solutions ip, ip* to the linear problems (3. IS). For continuous hierarchies, 
pseudo-differential analogues of the operators (6.16) and corresponding adjoint flows were studied in 

As an example, let u s show that taking ip, ip* to be the Baker- Akhiezer function tf>{l; z) and its dual, 
one reproduces eq. (|6.6| ). According to the prescription of Sect. 5.3, the Baker- Akhiezer function and its 
dual read 



l 7 a 

0S\l F Pz+l 



I , a 



-af\03\-l r P2-l 



I. a 



Pz=0 



(6.17) 



Pz=0 



where ( a is a normalization factor specified below. Here p z is the time variable corresponding to the 
flow Xqz left adjacent to I. Let z — > \\. In the limit p z coincides with u. Plugging ( |6. 17 ) with C, = X° z 2 
into (6.16), we do reproduce the operator in the r.h.s. of eq. (6.6). 



7 On the hierarchy of bilinear difference equations 

Integrable partial differential equations always can be included in an infinite hierarchy. Infinite number 
of commuting flows generates infinite number of evolution equations. 

The hierarchies of discrete integrable equations are less studied. First of all, it is not quite clear what 
are "higher discrete flows" on the space of pseudo-difference operators. An understanding of this matter 
is necessary if one is going to extend the Zakharov-Shabat formalism to the higher Hirota equations 
known in the literature |0|, fjof . 

There are two "complimentary" points of view on this matter. First, one might consider 3-term HBDE 



(4.20) as a counterpart of the whole infinite hierarchy itself. In this case it should be understood as an 



infinite set of equations (continuously numbered by labels A Q 's) for a function of infinitely many variables 
l a p associated with A Q A^. Second, one might expect that composite discrete flows are good candidates 
for true analogues of the higher continuous flows. This is justified by analysing the continuum limit. 
Indeed, to get a higher continuous flow as a limiting case, one should start from a composite discrete 
flow with specially adjusted labels. 

Our goal in this section is to show how these two approaches could be consistent with each other. 
A natural conjecture is that the iV-term "higher" Hirota equations for a function of N variables are 
consequences of the basic 3-term equation ( 4.20| ) treated as a hierarchy. This means that the 3-term 



equation is supposed to hold for each triple of the N variables with the corresponding A Q 's. To support 
this conjecture, the case of 4-term HBDE is considered in detail. Besides, an extension of the Zakharov- 
Shabat scheme to this case is suggested. 
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7.1 Higher equations of the hierarchy 



Higher Hirota difference equations known in the literature |5CJ are written for a function t(Ii, . . . , In) 
of N variables. They have the form 



1 z\ z\ 

1 Z 2 Zo 



„iV-2 



M-2 



nn 



T2T2 



1 Z N z 2 N 



N-2 
Z N T N T N 



= 0. 



(7.1) 



where zi are arbitrary constants and 

Tj = t(Ii, I2, ■ ■ ■ , h-i, li + 1, h+i, ■ ■ ■ , In) , 
n = r(h + 1, l 2 + 1, . • • , h-i + l,k, h+i + 1, . . . , In + 1) ■ 

In the more compact form they read 

N 

The constants Aj are subject to the only condition Ylj=i — 0- F° r ^ — 3 one gets usual Hirota's 
bilinear difference equation, where Zi are certain rational functions of A Q 's. Like in the 3-term case, the 
variables U will be identified with elementary discrete flows. 



(7.2) 



(7.3) 



The transformation 

tQi, exp 
sends eq. (f^) to the canonical form, 



N / n \ 

fc=i v=i,#fc / 



r(?i, . . .,x N ) 



N 



^TjTj =0, 



(7.4) 



(7.5) 



which does not contain any free parameters. 
An analogue of Hirota's original form, 



N 



^Aj-exp^Oj r-r = 0, 
is obtained from eq. (|T.S|) by the linear change of variables 



(7.6) 



1 N 

1 1 N - 2 ^ 



(7.7) 



i=i 



which generalizes (2.15). 



7.2 Zero curvature conditions for composite flows 

We are going to show that the zero curvature condition written for the composite discrete flows introduced 
in Sect. 3.1 lead to the "higher" bilinear equations of the form (7.1). The "higher" M-operators are 
obtained as products of the elementary ones. 

In this section we deal with the following graph of flows: 
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The reference flow is u. Other notation is clear from the picture. For simplicity, we consider left adjacent 
flows only. All what follows can be easily reformulated in terms of the right adjacent flows. 

According to the definition of composite flows (Sect. 3.1), we introduce a "higher" M-operator M qr 
generating evolution in the composit e flow labeled by the pair of vectors A0A2, A0A3 as the product of 
elementary M-operators of the form (4.15): 



Or 



(7.8 



or, in more details, 



M^(q, r) = M*(q, r + l)M r u (q, r) . (7.9) 

Due to the zero curvature condition for elementary flows M qr = M^ q . The compatibility of this composite 
flow with an elementary flow A0A4 reads 

M£(p, q + 1, r + l)Mf [p, q, r) = M qr (p + 1, q, r)M*(p, q, r) . (7.10) 



Clearly, t his z ero curvature condition is an immediate consequence of eq. ( |4.9| ) and definition (7^) pro- 
vided eq. (h9) holds for any pair of the flows from the triple p, q, r. 

A few words about the notation. For the notational simplicity, we put 



A' 



01 



A r 



A 1 



01 



A„ 



A 1 



111 



A. 



Since in this section we deal with too many variables, it is also convenient to change slightly the notation 
for the r-function writing the argument it in brackets: t u — > t(u). Other variables are indices, as before. 
This will remind us that u is the reference variable. Needless to say that unshifted variables will be 
skipped when possible. Throughout this section, the notation 



VM(p,q,r) = 



T P ' q ' T (u)T P+1 ' q ' r (u + 1) 
TP +1 '1' r (u)TP'1' r (u + 1) 



(7.11) 



is used. 



Now we are ready to elaborate eq. (7.10) explicitly. Eq. ( 7.10| ) reads 

(e a " - A P K (P) (p, g + 1, r + 1)) (e d « - A 9 K (9) (P, 9, r + 1)) (e a 
- (> _A,K (g) (P+l,g,r + l)) (e a "-A r K (r) (p+l,g,r)) (e 6 



\ r V^(p,q,r] 
\ p V^(p,q,r] 



(7.12) 



The coefficients in front of e 2a " give 

A P K (P) (P, q + 1, r + 1) + XgV^ip, q,r + l) + \ r V^%{p, q, r) 
= KVu+2 (p, 1, r) + X q V^ (p + 1, q, r + 1) + A r (p + 1, q, r) . 



(7.13) 



This is a direct corollary of 3-term HBDE. To see this, recall eq. ( 4.1q) from the proof of Theorem 4.1. 
In the present notation it reads 



A P K (P) (P, Q + 1, r + 1) + X q V^ (p, g, r + 1) = A p yj^ (p, g, r + 1) + A,K (?) (p + 1, g, r + 1) . (7.14) 
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Here r + 1 enters as a parameter. A similar equation can be written down for the pair p, r, where now 
q enters as a parameter: 



XpVj+M Q,r + 1) + X r v£> 2 (p, q, r) = X p V^%(jp, q, r) + KV^ip + 1, q, r) . 

Taking their sum, we obtain eq. ( 7.13| ). 

Comparing coefficients in front of e a ", we get: 

\\V^ (p, q, r + 1)K (P) (p, q + 1, r + 1) + A 9 A r y„ ( ;\ (p, g, r)v£\ (p, g, r + 1) 
+ A p A r (p, g + 1, r + l)V^ r + \ (p, g, r) 

= A P A 3 V^ (p + 1, g, r + 1)V^\ (p, g, r) + A 9 A r K (r) (P + 1,9, r)^ (p + 1, g, r + 1) 
+ A p A r l/ t +i (P- 9, r + (p, g, r) . 



(7.15) 



(7.16) 



In a similar way, it is easy to show that this equality also follows from HBDE. Of course this is a trivial 
result provided the zero curvature conditions ( |4.9| ) hold for both elementary M-operators in the product 
(©■ 

The fact which is not obvious from the very beginning is that eqs. (7.14), ( |7.15| ), (7.16) imply one of 
the higher Hirota equations (namely, (7.1) at N = 4). Proceeding like in the proof of Proposition 4.1, 
we rewrite (7.16) in the form 

t{u+1) ( T q+l ' r+1 {u + 2) TP +1 - r+1 {u + 2) 

A/iAr ~~; — 7~r~, r — AriA^ 



t(u + 2) \ q r r?+ 1 ^+ 1 (u + l) P r rP+ 1 - r + 1 (u+ 1) 



ApAg 



T P+hq+hr + l ( U )TP+ 1 > r + 1 (u + 1) T P+ 1 (u + 1) 



p+l,g+l,r+l 



(u+1) 



A„A r 



(u) 



r-r+l 



A„A 



(u) 



A„A 



r p+l,,+l,r+l( u ) ^ 1 r T P+l{u+\) P q T r + 1 ( U +l) 

^ +1 ' r+1 ( U )T(u + l)T r+1 (u + 2) 



(7.17) 



(cf. (4.16 



we get 



where 



p I V+ 1 (u + l)T-?+i.r+i( u +l) T (ii + 2), 
Multiplying both sides by (A 9 - A r )T« +1 < r+1 (u + 1)tP +1 (u+ 1) and using HBDE in the form 

{Xg ~ X r )T q+1 - r+1 { U )T(u + 1) = X q T r+1 { U )T q+1 (u + 1) - X r T q+1 (u)T r+1 (u + 1) , (7.18) 



T(U+ i)T p+1 ' q+1 ' r+1 (lL) 
Till + 2)TP+ 1 ' q + 1 ' r+1 (u + 1) 



yl p '«' r (u) = B p ' q ' r {u). 



(7.19) 



= (A, - X r )X q X r T q+1 - r+1 (u + 2)t p+1 (u + 1) 

T p+ ,q+ + 2) /x rJ)+ljr+1 ,, q+1>r+1 { 



An A, 



P 1 T p+l,q+l,r+l I 



r p+l,r+l 



> + l)) 



x , (u + 2) 

A P Ar T P+ 1 .9+ 1 .'-+ 1 (u) 



(X r T 



(U)T" +I > r+I (u + 1)), 



B p ' q ' r (u) = (A, - A I .)A g A r r p + 1 (u)r 9+1 ' r+1 (u + 1) 
ApA 9 ^^(A 9 rf+ 1 ( U + l)r q+1 (u + 2)) 



-+- Ap A r 



(«). 

r(u + 2) 



(A I .r p + 1 (u + l)r r+1 ( U + 2)). 
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The last two terms can be transformed further using eq. (7.18). The result is 



A™> r (u) = h p ' q ' r (u + 1) - X 2 p {X q - X r ) 
BP-i- r (u) = hP^ r (u)-X 2 p (X q -X r ) 



r g+l,r+l 



(u) 



r p+l,g+l,r+l 



(u+ 1)tP +1 (m + 2) 



T p+l,q+l,r+l( u j 

t +1 ^ r+1 {u)T{u + 1) T p +1 {u + 2) 



-(u + 2) 



where 



h p - q - r (u) = X q X r {X q - X r )TP +l {u)T q+l - r+ \u + l) + X r X p (X r - A p )r 9+1 (u)' 
+ A P A,(A P - A 9 K +1 ( u )tP +1 ^+> + l) . 



.p+l,r+l 



(U + 1) 



Finally, eq. ( |7.19|) turns into 



-(u + l)T p+1 ' q+1 ' T+1 (u) 



hP^ r {u) 



t(u + 2)TP+ 1 'i+ 1 ' r+1 (u + 1) hP>i> r (u + 1) ' 

which leads to 

h p ' q ' r (u) + (A p - A,)(A, - A r )(A r - X p )t(u + l)r p +^ q +W («) = . 
Eq. ( 7.22| ) has the form (7.1) at N = 4. This completes the calculation. 



(7.20) 



(7.21) 



(7.22) 



Af-operators for arbitrary composite flows can be defined as a straightforward generalization of (7.8): 

(7.23) 



/V 



N 



MP* 



«p(EMlI( e ~* lM ^ 



Note that the order of the operators in the product is not essential since the operators e~ dp i commute 
due to the zero curvature condition (LE). For simplicity it is assumed that all the flows pi are left adjacent 
to the reference flow u. Operators ( 7.23 ) generates discrete analogues of higher flows of the KP hierarchy 
(see the corresponding graph of flows in Sect. 3.2). There is now straightforward to write down similar 
operators for right adjacent flows which would generate higher flows of the discrete 2DTL hierarchy. 



Conjecture. All the higher HBDE (7.1) follow from the compatibility of the composite flows gen- 
erated by the M-operators ( 7.23| ) and elementary discrete flows. 

The calculation given above shows that the conjecture is true for the 4-term bilinear equation. 
Unfortunately, we are not aware of any other proof than this sophisticated calculation that is hard to 
perform in the general case. The conjecture claims that all the higher bilinear equations are corollaries 
of eq. (4.2C) considered as a hierarchy, i.e., valid for all triples of adjacent flows. 



8 Reductions of Hirota's equations 

The hierarchy of discrete Hirota's equations admits several important reductions. Reduction means 
imposing a constraint compatible with the hierarchy, so that the number of independent variables gets 
reduced. In this way one is able to construct discrete analogues of KdV, SG and other interesting 
equations. 

The simplest way to impose a constraint is to require the r-function to be stationary with respect to 
a particular flow s (possibly, up to a "gauge" transformation ( |2.2| )). Non-trivial examples emerge when 
the stationary flow is composite. As for the commutation representation, there are two options. 

First, the stationary flow can be reference. Then M-operators become free of differentiation because 
the symbol d s commutes with all the coefficients. In other words, d s can be considered as a c-number 
and identified with a spectral parameter. This is a natural origin of Ai-operators depending on (rational) 
spectral parameter. 
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Alternatively, one may take as the reference flow any flow y other than s. Then any M-operator 
y[ (/) generating a flow / contains operators d y . Since coefficients of the operator do not depend 

on s, the compatibility condition for the flows s and / acquires the Lax- type form 

Af (s) (/ + l)M (/) = M {f H&\f) , 

where JVf" ' plays the role of Lax operator. Unlike the Zakharov-Shabat scheme, where each zero curva- 
ture condition involves two different time flows (apart from the reference flow), Lax equations are written 
for each flow separately. The Lax equation represents the time flow as a similarity (i.e., isospectral) trans- 
formation of the Lax operator. It is natural to call the L-operator of the reduced hierarchy. This 
is a natural origin of L-operators which are difference (or differential) rather than pseudo-difference (or 
pseudo-differential) operators?]. In the light of this general scheme, we pass to considering examples. 



8.1 KdV-like reductions 



0). Discrete d'Alembert equation (a trivial example). Let u be an elementary discrete flow dealt with in 
Sects. 3,4. The stationarity condition with respect to this flow, r(u + 1) = t(u), immediately leads (see 



flh20| )) to the relation 

where I, m are any other two elementary flows. This is the discrete 2D d'Alembert equation written in 
"light cone" coordinates. The general solution is T l ' m = x+(l)X-( m ) with arbitrary x±\ this is just the 



allowed "gauge" freedom (2^8) of the r-function. So we see that this reduction is too strong - it respects 
trivial solutions only. In this sense it is incompatible with the hierarchy. To get nontrivial examples, one 
should either impose stationarity conditions with respect to the composite ("higher") flows or impose 
periodic conditions in u with period N > 1 (e.g. t(u + 2) — t(u)). 

1). Discrete KdV equation. Consider the graph of flows 




and put A° 



01 _ \ \01 
A <3> A 3 



X p for brevity. In this notation, eq. ( 4.2C| ) takes the form 

A g T„ T u+l — A p T u T u+1 -(- (A p - A q )T u T u+1 — U . 

To get the discrete KdV reduction, we impose the constraint 



U U 3 



(8.3) 



i.e., the r-function is stationary with respect to the composite flow labeled by the pair of vectors A0A2, 
A0A3. This condition converts the 3-dimensional equation (8.2) into the following 2-dimensional one: 

Krl^Tlll - \prl- l T p u X\ + (K ~ K)t P u< +1 - . (8.4) 

This is the discrete KdV equation in the bilinear form Q , [Q . The discrete KdV equation is also known 
in the form |3(J 

V(u,p) - V(u- l,p- 1) = k(V _1 (u,p- 1) - V-\u- l,p)). (8.5) 

6 The Lax operator for the hierarchy without any reduction is in general an infinite series in negative powers of first 
order difference operators. The theory based on the Lax representation with Lax operators of this kind is not considered 
here. 
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The equivalence to eq. (8.4) follows from the identification 



V(u,p) = 



~u~u+l 
_P+1_P 



A,, 



Let us turn to M and L operators. Following the history of the KdV equation, we begin with 
difference operators with scalar coefficients. Let u be the ref erence variable. Then the (composite) 
stationary flow is generated by the M-operator of the type (^) : 

^ 'u ' U + l , x 'u + l 'u- 

A P p p-1 + A 9 p-1 p 



MJ9 = M£(p, q + l)M«(p, q) = = e 29 " 



'« 'u+l 



ApAg . 



6) 



u+l 'u+2 , 



We call this 2-nd order difference operator the Lax operator of the discrete KdV equation. The spectral 
problem L^ KdV) ip = Ei/j is a discrete analogue of the stationary Schrodinger equatio n that is an auxiliary 
linear problem for KdV. The p-evolution generated by the M-operator of the type ( 4.1 5| ) is isospectral: 
L (KdV)( p + ^ _ MPL( Kdv \p){MP)-\ 

If the reference flow coincides with the stationary one, we get 2x2 matrix M-operators depending on 
a spectral parameter z. The spectral parameter is an eigenvalue of the shift operator along the stationary 
flow acting to the ^-function: exp(<9 p + d q )ip = z 2 ip. Consider the vector 



V xl J 



r P+l 



V 



.p+i 



J 



Repeating the argument of Sect. 5.5, we obtain the following linear problems for shifts in u and p 
respectively: 



( €+i \ 



V xl+i ) 



T P T P+1 
'u 'u+l 

~P+ 1 -P 

U 



v+1 



V X P u +1 J 



V 



v+1 



'u+l 



r P+2 
' u 

p+1 



1 1% \ 



J 



i.7) 



jp+i 



V Xl J 



Xn-X 



P J 



\ xl J 



Compatibility of these linear problems yields the discrete matrix Zakharov-Shabat equations with the 
spectral parameter z. 



Equations for the first component tp% read: 



€+i 



1 



Xp — X q 



_P+1 P-1 
,2 'u-1 'u+l 
K q _P _P 
'u-1 'u+l 



-X 



p-1 p+1 ' 
2 'u-1 'u+l 



P -P 



V u = (z 2 - A P A 9 )^_! , 



'u-1 'u+l 



^ +1 + (A P -A g ) p ^' p+1 ^ u = z^. 



p\2 



(8.9) 



(8.10) 



Eq. ( |8.9| ) coincides with the spectral problem for the Lax operator (8.6) provided the r- function obeys 
the bilinear relation (|Q|). 

2). Discrete Boussinesq equation J7^ /. In this example the T-function T p,q satisfies T^\ ,q+1 — tP ,<! . 
Eq. ( pT2] ) reduces to 

XqT P+l,q T p-l,q _ XpT P,q+l T P, q -l + _ X ^ T P+l,q+l T p-l,q-l = Q > (g U ) 

where u enters implicitly as a parameter. Different kinds of Lax and Zakharov-Shabat representations 
for this equation can be written down straightforwardly. 

In a similar way, it is possible to define A„-type reductions. In this case the T-function obeys the 
condition exp ( X}«=i dp a ) T = T - This becomes an actual reduction for higher Hirota equations with the 
number of variables greater than n. 
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8.2 Discrete time ID Toda chain and its relatives 



This group of examples includes the discrete time ID Toda chain (1DTC), discrete AKNS system (in 
particular, discrete nonlinear Schrodinger equation), discrete time relativistic Toda chain and discrete 
Heisenberg ferromagnet (HF). These models differ by choosing dependent and independent variables 
while the type of reduction is essentially the same. 

Let the graph of flows and the notation be the same as in Sect. 4.3. Now the r-function is required 
to be stationary with respect to the composite flow labeled by the pair of vectors A0A3, A1A2, i.e., 



T l+i,l _ T u 



.12) 



The stationary flow is generated by the "composite" M-operator e dl M l n e dl M l n with the reference flow 
n. This operator should be identified with the Lax operator of the discrete time 1DTC: 



L (TC) = 



'n'n+l 
'n T n+1 



T l + 1 T l 
1 tJ-tH-I 



r ri-l T n+l - 



.13) 



This is a 2-nd order difference operator in n with scalar coeffi cients. The spectral problem L^ TC ^rp = 
Ei\) is a ID discrete stationary Schrodinger equation (cf. ( 5.43 )). The Z-dynamics does not change the 
spectrum of L^ TC \ 



Replacing e c 



by 2x2 matrices depending on the spectral parameter 



we get from ( 4.23 ) an L-M pair for the discrete time 1DTC realized 



/ t 1 t 1+1 



II- 



'n+1 



T n+1 \ 
-llV—j- 
T 



'n+1 



.14) 



M^\z) 



V 



T n+1 

-t^-iZT 



2 'n+l'n-1 
-Z — (1- 



In T 1 



(8.15) 



Correspondingly, the 2D discrete Schrodinger equation (5.44) becomes a ID spectral problem in the 
variable I, 



1 



TlT, 



n'n+l ,1 



i.l 



V- 



r l-l T l 

'n 'n+1 .1,1-1 



' " ' n- I \ T n T n+l 



^n 



(8.16) 



with the non-standard dependence on the spectral parameter. Another useful form of this equation is 



1+1 



'n 'n+ljl-l 



■n-\-i it — 1 It/ \ — 1 n n+1 i / 



T n T n+l 



T n T. 



.17) 



n+1 



where t/c = z 



The equation of motion in the bilinear form is the 2-dimensional reduction of eq. (2.11): 

= iy - m)(t„) 2 , (8.18) 



v ' n ' n 



which after the linear change I — > rn = I + n, r l n — > r„(m) coincides with Hirota's original form of this 
equation ||: 

(1 +g- 2 )T n (m + l)r„(m - 1) - t„ + i(to)t„_ 1 (to) = g- 2 (T n (m)) 2 , (8.19) 



where 



g- z = --i- 

A 4 
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In terms of the new dependent variable 



0„(m)=log^±iM (8.20) 
T n (m) 



eq. ( B.19 ) acquires the form 



U i i u oa f \\ 1 +5 2 exp(0„+i(m) - <f> n (m)) 

exp (0„(m + 1) + - 1) - 2*,(m)) = t +g2fflcp(0n(Tn) _^_ l(m)) ( 8 - 21 ) 



that is the discrete time 1DTC equation studied by Yu.Suris |31 



The continuum limit in m is straightforward: set m — > ro/e, g 2 = — e 2 , then (j> n (m± 1) — > </>„ it^J, + 
ie 2 0" as e — > 0. Expanding eq. (3.21) in e, we get the familiar 1DTC equation: 

0" = e ^-*"-i - e *"+ 1 -*" . (8.22) 



It is interesting to note that eq. ( 3.21 ) has another continuum limit yielding the sine-Gordon equation. 
To see this, let us redefine the field <j) before taking the limit: 

ct> n {m)=i{-l) n - m cp n {m), (8.23) 

so the equation reads 

/ / \ i \ i i\ i \\ 1 + 9 2 ex P(*</'n+l( TO ) + *V?n( m )) so n/\ 

exp [iip n+1 (m) + itp n -i(m) - iip n {m + 1) - i<p„[m - 1)) = — = -. — : r : — r — -. (8.24) 

1+g ' exp{-i(p n -i(m) - tip n \m)) 

It is the field ip n (m) that is now assumed to have a smooth continuum limit in n, m. Setting m — > m/e, 
n — > n/e, g~ 2 = — ie 2 , e — > and expanding in e as before, we get the SG equation 

(^-^V = 2sin(2<^). (8.25) 

(The same limit for the field in eq. ( |3.2l| ) would give the d'Alembert equation (9 2 — d^ip = 0.) A 
discrete analogue of the SG equation of a different kind is given below in Sect. 8.3. 

Remark 8.1 The discrete KdV reduction discussed in Sect. 8.1 is formally a particular case of the 
present reduction when A and Ai merge. However, it is more convenient to consider them separately. 
Note also that there exists a continuum limit leading to the KdV equation from eq. ( |8.18 ), too. In this 
sense the discrete time 1DTC is sometimes considered as a discretization of the KdV equation. 

The discrete version of the decoupled nonlinear Schrodinger equation (called also the AKNS system) 
is essentially the same reduction, i.e., the stationary flow is the same. The difference is in picking out 
other independent flows to be involved in the equations. Specifically instead of the n-flow one may 
consider any other elementary discrete flow p left adjacent to I. 

We would like to show how to derive the discrete AKNS system directly from the bilinear Hirota 
equations. The basic bilinear equations in question are 

\ P rf +1 T^+(» X^rfr^ 1 = rift 1 !*** , (8.26) 

- V£M:F +1 = {u- \ p )tZ + ^1? . (8.27) 

The first one is eq. (4.22) written for the triple (p, n, I). The second one is obtained from eq. ( 4.20) ) 
written for the triple (l,p,n) by taking into account the invariance r '+ 1 - / + 1 = r M It is easy to see that 
in terms of the quantities 

^pj _ -P-J 



qp,i = _M£ R p,i = ( 8 .28) 

i n ' n 

(n is fixed) these equations can be rewritten as follows: 

(v - A p )(A p - fi)QP J+1 = (Ap - it Q pJ+1 R p+lJ ){v Q p > f - X p Q p+lJ+1 ) , (8.29) 
(f - A P )(A P - f i)R p+1J = (A p - liQ pJ+1 R p+1 ' T )(vR p+1 ' 1+1 - X p R pI ) . (8.30) 
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This system is equivalent to the discrete AKNS system given in ref . |13[ . 

There exists another choice of dependent variables which converts eqs. QS.18 ), (8.26), ( [B.27 ) into a 



discrete analogue of the relativistic Toda chain (RTC) J52|. Passing again to the variable m = l + n, like 
in eq. ( |3.20| ), we put 



x m {p) = - log 



rg(m + l) 
T«(m) 



.31) 



The claim is that the equation of motion for x m {p) has the form 

(l-aexp(x m+1 (p)-x m (j))))(l-(3exp(x m (p)-x m -i(p - l)))(l-^exp(x m (p)-x m {p + 1)) 
(l-aex.p(x m (p)-x m -i(p)))(l-(3exp(x m+ i(p + l)-x m (p)))(l-jexp(x m {p - l)-x m (p))) 

where 



= 1 . 



An 



32) 
33) 



M - A* A p (/i - v) ' 

which is the discrete time RTC equation suggested in ref. |33j. 

Let us outline the way to derive eq. ( B.32 )^[ The basic bilinear relations ( B.18 ), ( 8.26 ), ( S.27 ) read 

(u ^(m)) 2 - vrl{m + l)r^(m 1) = -^(m^H , (8.34) 

ifi \)i*{m + 1K +I (m) + \ p rl +l {m + l)<(m) = K +1 (m + lK^(m) , (8.35) 
< +1 (m + 1)< +I (m - 1) - \ p rl{m)rlX\{ m ) = (y \)r^ x (m)^ +1 (m) , (8.36) 

respectively. It is straightforward to show that the following two bilinear relations are direct corollaries 
of the basic ones: 



(" - ^X+iMrr^m + 1) - ( M - \ P Ktl(mK(m + l) = (u- ^)r p n+1 (m + l)r£+V) , 



.37) 



\ P { V - nK(m)T? +1 (m) + K/i - A P )r£(m + 1K +I (m - 1) = - A p )< +1 (m)^ (m) . (8.38) 



r p+i 



p+i. 



Eq. (8.37) is obtained by eliminating T^(m) from ( 8.35j ), ( 8.36| ) (i.e. by dividing them by r^ +1 (m + 1) 



T n+i( TO ) respectively and taking their sum) and making use of eq. ([B~3ll). Eq. 
similar way by eliminating r^ +1 (m) from ( |S.35 ), ( B.36| ) and making use of eq. 
easily follows from (|3§), fl8.35| ) and ( |8.38| ). 



.38) is obtained in a 
157]). Now eq. ( |8~3^ ) 



Remark 8.2 The bilinearization of the usual (continuous time) RTC was suggested in ref. j34|. The 
equivalence of the RTC and the "semi-discretized" AKNS system (with discrete "space" and continuous 
time variables) was recently proved in ref. |35[ ]. 

We conclude this subsection with a few words about the discrete Heisenberg ferromagnet [|l2j . This 
equation fits the scheme in the following way. The reduction is the same. In addition to the flow p 
from the previous example, one should introduce yet another elementary discrete flow q left adjacent 
to I. The r-function now depends on four variables: r = T^(p,q). Fix n, I and consider the following 
four functions of p, q: r^ +1 (p, g), r^ 1 (p,q), r^t* (p, q), t„1i (?>)<?)• It can be shown (using, in particular, 
one of the higher Hirota equations in 4 variables) that certain combinations of these functions satisfy a 
system of nonlinear difference equations in the variables p, q. This system is equivalent to the discrete 
HF model discussed in detail in r cf. @, where it was treated in a slightly different way as a part of the 
reduced 2-component 2DTL hierarchy. As in the case of the discrete AKNS system, the aforementioned 
embedding into the 1-componcnt 2DTL hierarchy leads to equivalent equations of motion. We omit the 
details. 



8.3 Periodic reductions 

Periodic reductions of the continuous 2DTL hierarchy give rise to a number of very important equations. 
For example, the 2-periodic reduction t„ + 2 = r„ contains the sine-Gordon (SG) equation. The same 
periodic constraint can be imposed in the discretized set-up, thus providing us with a discrete analogue 
of the SG equation. 

7 The idea of this derivation is due to S.Kharchev (unpublished). 
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We remark that periodic reductions can be treated on equal footing with stationary reductions. 
Indeed, the flow p — » p + 2 is, formally, a degenerate case of a composite flow when the corresponding 
labels pairwise merge on the complex plane. The periodicity r p+2 = t p means stationarity with respect 
to this degenerate "composite" flow. However, this point of view does not seem to be useful in practice. 
Usually it is more convenient to treat periodic reductions separately. 



Let us consider the 2DTL-like form of HBDE (2.11) with the constraint 



<i 2 = r^. (8.39) 
It becomes the following system of two-dimensional equations: 

i.i+i i+ij i \ i,i l+i.l+i l.l+i l+i,l 
vt q t - (y - (j,)t ' t = ht{ t x ■ , 

1,1+1 1+1,1 i \ 1,1 1+1,1+1 1,1+1 1 + 1,1 (Q / , n ^ 

vt{ rj - (y - h)t{ Tj = /ir r . (8.40) 
The SG field & l,L on the square lattice (I, I) is given by 

^ J =l l0 ^r ( 8 - 41 ) 

Rearranging eqs. ( 8.40| ), one gets a closed equation for 

v sinh($^ r + s'+W - $^+! - s'+i.i) = /x sinh($' >' + + + , (8.42) 

which is known as discrete SG equation Q written in the light cone coordinates. 
Let us mention another useful form of the discrete SG equation |36|, j37j. Set 

S l,l = cxp ( _ 2 $i+U _ 2 $M+!) = T ^^ r j'^| , (8.43) 

T o +1 " T o +1 
i,i i+i,T+i 

S 1 ' 1 = exp ( - 2$^ - 2$' +lJ+1 ) = T ] T T ] +1 , (8.44) 

then, evidently, 

gl,T+l jgl+i,l _ oMe/+i,'+i 
On the other hand, the discrete SG equation implies that 



/i - vS 1 ' 1 
fiS 1 ' 1 — v 



(8.45) 



so eq. 



S.42) converts into 



lWc(+1 , w _(/i-^W)^-^Hl) 



S ' 5 ' "(^M+i-^C/i^-^) (8 - 46) 
(cf. ©). 

We now turn to the zero curvature representation. Let I be the reference flow. The shift n — > n + 2 
is generated by the scalar L-operator 

(1,1 1+1,1 1+1,1 l+2,l\ 

However, this representation is not convenient for describing evolution in I. 

The matrix L-M pair with spectral parameter is more appropriate here. To derive it, one should 
take the "stationary" flow 2n to be the reference one and repeat the arguments given in Sect. 5.5 with 
necessary modifications. The operator e 29 " should be substituted by a spectral parameter z 2 . Omitting 
the details, we present the result. 
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The auxiliary linear problems read 



V x l+lj ) 
V x lJ+1 J 



1,1 1+1,1 

T o T i ' 
1+1,1 i,i 

T T l 

1+1,1 

Z T, 



1,1 
v M 



v 1+1,1 



,1+1 \ 



V 



z -M+i 



i,i 1,1+1 
T a T i 
1,1+1 i,i i 



( ^ \ 



\ x lJ J 



\ x 1 ' 1 J 



(8.49) 



'i 'o 

which is akin to (8.7). Denoting the matrices in the r.h.s. of eqs. ( B.48] ), ( 8. 49] ) by AfW, respec- 
tively we can write the compatibility condition 



M M (l,l + l)M (-) (/,/) = M(- ) (l + l,l)M {+) (l,l) 



(8.50) 



which yields the discrete SG equation. 

iV-periodic reductions (t„ + tv = ^Vi) can be treated in a similar way. They correspond to ./V-periodic 
Toda lattices in discrete time. It is also possible to impose periodicity with respect to any one of the 
composite flows. In the remaining part of this subsection we briefly comment on an important class of 
such reductions which are discrete analogues of Intermediate Long Wave (ILW) equations. 

The universal form of reductions from the 2DTL to the family of continuous ILW equations is most 
transparently written in terms of the r-function of the 2DTL hierarchy. The reduction to the ILWfe 
equation reads |38| 

T n +k(ti + h,t-i, ■ ■ ■ ; ti, ?2, • • •) = r n (ti,t2, ■ ■ ■ ; h, h, ■ ■ •) , (8.51) 

where h is a fixed parameter. This parameter interpolates between the fc-periodic reduction (h = 0) 
and the Benjamin-Ono equation (h —>■ oo). In words, the r-function should not depend on a particular 
combination of n and t\. This suggests a discretization of the ILWfe equation. According to our general 
rules of discretization, one should substitute t\ by an elementary discrete flow p. Then it is natural to 
substitute eq. ( B.51 ) by 

<Xl = Tl, (8.52) 

where I and k are integer parameters. The particular cases are the discrete KdV equation (I = k = 1) 
and the fc-periodic reduction (2 = 0). In the continuum limit we get the continuous ILWfe equation. 



8.4 Discrete Liouville equation 

The discrete Liouville equation (DLE) and its ^-generalizations J39| (discrete time 2DTL with open 
boundaries) form a very important special class of discrete integrable systems which in general does not 
fit the reduction scheme discussed in this section. We include it here for the only reason that the DLE is, 
formally, a degenerate case of the discrete SG equation. The relationship between these two integrable 
systems deserves further study. 

The DLE can be obtained fro m the discrete SG equation as a result of a certain scaling limit. 
Let us rescale S 1,1 — > ^iS 1,1 in eq. (8.46). Clearly, this rescaling means a constant shift of the field: 



<£>M <j)M _ i log /i. Then, taking the limit ji — > in eq. (|S.46) (keeping shifts in I alive!), one arrives at 



4 

the DLE 



S U S 1+1,1+1 = {v -i _ s U+i ){y -l _ gl+1,1) ( (8-53) 

where, formally, 

S l £ l = hm (m" 1 ^' 1 ) . (8.54) 



3G 



Setting 



S l £ l = cxp ( - ~ 2<t> l i I+1 ) , (8.55) 



we get, in place of eq. (8.42), the DLE written in terms of the discrete Liouville field 



2v sinh(<i>^ + <S> l + u+1 - - = exp + $ l + u+1 + 4>^ +1 + , (8.56) 



or, in a simpler form, 



exp ( - 2<S> l + U - 2<S> l i +1 ) - exp ( - 2$'/ - 2$ i + a+1 ) - v~ x . (8.57) 

In the continuum limit one should put I — > v 1 ^ 2 x+, I — > v 1 / 2 X-, S£ — > exp ( — 4$(x+, £_)). 
Expanding in z/ -1 — ► 0, we get the continuous Liouville equation 

2d x+ d x _<£(x+,x-) =e i *(*+> x -K (8.58) 
The bilinear form of eq. ( 8.53| ) is available via the substitution 



j,T_ T 1 (l + lj)T 1 (hT+l) 
T°(l + lJ)T 2 {l,l + l) 



s« = 'iT, (8.59) 



after which the DLE becomes equivalent (up to a " gauge freedom" , see below) to the bilinear relation 

T a {l + 1, l)T a (l, 1 + 1) - T a (/, l)T a {l + 1,1+ 1) = v- x T a - x {l + 1, l)T a+1 (lJ+ 1) (8.60) 
with the condition 

T a (2, 1) = unless a = 0, 1, 2 . (8.61) 

This condition implies the discrete d'Alembert equation (^l]) for T° and T 2 , so they have to have a 
factorized form T°(l, I) — x°(0x°(0> T 2 (l, I) — x 2 (0x 2 (0 with arbitrary and independent functions x°' 2 > 
X 0,2 . This is just the aforementioned gauge freedom. 

The striking similarity between eqs. ( |B.4g ) and (8.56) is transparent after the replacing T a (l, I) — > r l a ' 1 . 



Furthermore, taking into account the periodicity r 2 ' = r ' , they become formally identical. (Equiva- 
lently, using the gauge freedom, one can set T 2 (l,T) — T°(l,T) in eq. (8.56).) It would be interesting to 
link them directly on the level of solutions. 
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Appendix: bilinear difference equations from continuous hierar- 
chies 

In the Appendix we give an alternative point of view to the difference Hirota equations. It relies on 



the famous Miwa transformation (3.3) which so far was obscure in our exposition. Given a continuous 
integrable hierarchy (such as KP or 2DTL), this relation can be used as a definition of the elementary 
discrete flows. This definition leads to the same discrete flows as in Sect. 4.2. This approach has as many 
advantages as disadvantages. The main advantage is much more direct and instructive connection with 
the grassmannian approach to continuous hierarchies and their r- functions. The main disadvantage is 
misleadingly less invariant formulation which is inconvenient in some cases. 
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The Miwa transformation 



Let r(ii, t2, £3, . . .) = r(t) be the r-function of the continuous KP hierarchy. It is a function of infinite 
number of "times" U and satisfies infinitely many bilinear equations. The r-function solves all equations 
of the hierarchy simultaneously. 

In general the r-function can be represented as an infinite dimensional determinant |p^| - |l7| . It 
turns out that there exists a choice of independent variables such that the determinant reduces to a 
finite dimensional one. This choice is provided by the Miwa transformation [ pi] : 

**=*i 0) -^E^a*' k = l,2,... (Al) 

Here, the summation runs over a finite set /, are "background values" of the times, /i Q are arbitrary 
complex numbers (called Miwa's variables) and p a are integers (sometimes called multiplicities of \x a 's) . 

Remark The Miwa transformation plays an important role in revealing integrable structures of 
matrix models of 2D gravity. In particular, the easiest proof of the fact that the partition functions 
of the Kontsevich model [|t| and its generalizations are r-functions of the KP hierarchy relies on 
Miwa's transformation. 



In what follows we use the short hand notation fl3.q) . 
Important fact: The r-function of the KP hierarchy obeys the identity 

r + e (k- 1 ] - k- 1 ])) = n ^ r r )n °d a 7 g) ^ 1 <su k ^^ (A2) 

where 

Kiv a) (A3) 
Here N > 1 and /i Q , v a are arbitrary complex numbers. A useful particular case of this formula is 

r { t (o) - fx 1 ]) = dct ;r f MM) , (A4) 

V tt=l / ria</3(Ma -Hp) 

where 

<PM = 7 rry Hm v a — 1 — -j-iif ( v , M ) . (A5) 

(m — l)! i/^oo yi/ m 1 



When one translates the KP theory into the language of free fermions [[L5| formula (A2) becomes nothing 
else than Wick's theorem, K(y, fj,) being the fermionic propagator on a Riemann surface. 



Instead of treating eq. (A2) as an identity, one may follow another way. Given a function K(v^i) 
with a simple pole at v — fj,, this equation can be used as a definition of the l.h.s. This simply means 
that we disregard the dependence on background times vP assuming they are fixed. The r-function in 
Miwa's variables satisfy certain bilinear relations to which formula (A2) gives a solution in the form of 
a finite dimensional determinant. 

In the case of the 2DTL hierarchy the Miwa transformation goes in a similar way. The r-function 
T n (ti,t2, ■ ■ ■ ; ii, fa, . . .) = T n (t;t) depends on the discrete time n and two infinite sets of continuous times 
ti and ii. We set 



^ = 4 0) -^E>^' k = l,2,... (A6) 
where p, a is an independent set of Miwa's variables with multiplicities p a . 



The following analogue of eq. ( A2) holds 

/ N N \ 



Tn-N [ t<» m + J>j ) = 'f ( /' (0))n f_ 1 ^"\ det ; Jn^fifi), 



Ct = l 



ce=l 



where 



Jn(M,M) 



Note that in this case the function J n (/j,,ji) does not necessarily have first order pole at fi = p,. 



(A7) 



(A8) 



Discrete flows 

Discrete equations for the r-function listed in Sect. 2 are obtained if one fixes Miwa's variables fj, a and 
consider flows in the multiplicities p a . We give a few examples. 



Example 1. Set 



r Pl,P2,P3 — 



Tit 



(0) 



(A9) 



and consider = f(°) — X^q^i-Pq^q 1 ] as a new "background". According to eq. (A4), 

r^ +1 =^( m ), 



Pl + l,p 2 + l _ 



^i(Mi) Vi(f^) 



Ml - M2 

with some functions (f>2- Combining the zero determinant with two identical lines, 



(A10) 





^i(Mi) 


^i(Ma) 


^l(M3) 


= 


0i (Mi) 


^i(Ma) 


^l(M3) 




<^2(Ml) 


£2^2) 


^2(M3) 



(All) 



and expanding it in the first line, we get an equation of the form (2.7). Since its coefficients do not 
depend on the chosen background, the equation holds for all values of Pi,P2,P3- 

Example 2. Repeating the previous argument for 

= r (V» +po([^ 1 ] K 1 ]) + X>«([" _1 ] - [Ma 1 ])) 
and making use of eq. (A2), we get eq. ( [2.10 ). 



Example 3. In eq. (A7), let N = 2, fix = fi, fix = fi, /12 — > 00, /i 2 — » 0: 

r„_ 2 (t(°)-[M- 1 ];? 0) + [M])T„(t (0) ;t (0) ) = 



r n _!(t(°) - [ M - x ]; <W + [mD r„_i(*(°) - [m- 1 ]; t (0) ) 



Denoting 

we get the equation 



'n ~~ 



(tW-iD*- 1 ];^-^]), 



,/,J+M+l,i _ 1,1 1+1,1+1 _ f -/ X ij+l i+l,i 

'n 'n — UVWn+1 1 n-1 



n n 



(A12) 

(A13) 
(A14) 



Example 4- Example 1 can be generalized in the following way. Consider a N x TV-matrix with the 
lines (pi(iii), cpi(fii), (p2((ii), <p${pi), tpN-i(fXi), i = 1, 2, . . . , N, so that the first two lines coincide 



39 



and determinant of this matrix is zero. Then, expanding in the first row, like in Example 1, we get the 
"higher" bilinear difference equation of the form (^J). 

Example 5. At last, we show how to derive HBDE in the KP-like form from eq. (A4) in a direct wayfj. 
When two or more jU a 's coincide, both the numerator and denominator in the r.h.s. of eq. (A4) equal 
zero. Resolving the indeterminacy, we have 



T t 



(0) 



N \ 
J^Palpa 1 ]) = 
a—1 / lla< 



det (Mf >) 



l<i,j<M, 



and all fi a 's are now distinct. Here JV = Yla=iPa> ^ij^ * s * nc ^ x A/"-matrix having the rows 

<Pi(t*i), Vi(Mi)> </>"(Mi): • • • > ¥ J i Pl_ (Mi)> 
^»(M2), Vi(M2), Vi'(M2), • • • , ^ P2_1) (M2), 



(A15) 



We need the well known Jacobi identity for determinants: 

D[ii\j 1 }D[i 2 \j 2 ] - D[ii\j 2 }D[i 2 \ji\ = D[i 1 ,i 2 \j 1 , j 2 ]D , 



Ki<Af. 



h < *2, ji < h , 



(A16) 



(A17) 



where D is determinant of a square matrix and D[ii, i 2 \ji,j 2 ] denotes minors of this matrix with ii^-th 
rows and ji^-th columns removed. Applying this identity to the matrix M^j in ( A15) for i\ = Af — 1, 
H = A/\ ji = Ysl^iPa, h = h + ELa+iP«> 1 < a < 6 < A, we get, in the short hand notation, 



(fJ-a ~ Hb)TT 



Pa-l,Pb-l _ _Pb-l £P, 



= r" r — r t 
V(JV-i) 



M 



(A18) 

for 1< i < AT -2, 



where f is defined by the same formula ( A15 ) with the matrix M q 

M {N ~ 1] - M (A ° 
Af—l,j ~ M Af,j ■ 

Let /i c be a third Miwa variable (different from /x , /xj) with the multiplicity p c not shown explicitly 
in eq. ( Al8| ). Multiplying this equality by t Pc ~ 1 /t and then writing down a couple of similar equations 
obtained by cyclic permutations of the indices a, 6, c, we see that the sum of these three equations 
coincides with eq. fl2.7p . 

Remark The discrete flows discussed here coincide with those introduced in the main body of the 
paper if one fixes the following choice of the labels Ao and Ax: Ao = oo, Ai = 0. (To remove a label to 
infinity, one should use a different normalization.) 



Continuum limit 



As it is clear from eq. (Al), inverse Miwa's variables /i" 1 play the role of lattice spacings for the discrete 
flows. So, to perform the limit to continuous equations, it is necessary for fj, a to tend to infinity with a 
simultaneous rescaling of p a 's. 

Here is the typical example (the KP hierarchy). In this example we follow ref. (3^]. Introduce three 
(a priori independent) lattice spacings £{ — \ij , i = 1, 2, 3, and rescale pi — * pt/si. It is then convenient 
to rewrite the KP-like form ( fj77| ) of HBDE in terms of Hirota's D-operator (1.2): 

(ei{e 2 - e 3 )e~ (£l/2)D ^ +(£2/2)D ^ +(£3/2)D ^ + 

+ £2 ( £3 _ £l ) e (ei/2)A>i-(^/2)A> 2 +(W2)A, 3 + 



+ e 3 ( ei _ £ ^ e (ei/2)D P1 +(e2/2)D P2 -(e 3 /2)D P3 ^ 
3 this argument is taken from ref. J4l|], 



T • T = 0. 



(A19) 
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This equation serves as a "generating function" for a part of the continuous KP hierarchy. To see this, 
we express operators D Pi through Hirota's derivatives with respect to the continuous flows t k , 



oo 



i= 1,2,3, 



k=l 



in accordance with (Al). Substituting this into eq. (A19) and expanding it in a power series in £j, we 
have: 

ei(e 2 -e 3 ) £ e {el E \ Vj {] i b)Vu{-\b)Vi{-\b)+ 

j,k,l=0 

oo 111 

j,k,l=0 

°° 111 

es(ei-e 2 ) ]T 4 £ *e l 3 V j (--D)V k (--D)V l (-D) ] r • r = 0, 

j,k,l=0 



(A20) 



where D = (D tl , D t2 /2, . . . , D tk /k, . . .) and Vj(t) are Schur polynomials defined by 

(oo \ oo 

k=l I m=0 

Extracting the coefficients in front of £^£263, we obtain the infinite set of bilinear equations, 

TVi(l-D) Vj-i(-iD) Vj(r\D) 

v k -x{\b) v h -x{-\b) v k {-\b) t-t = o, 
Vi-i(±b) Vi-^-^b) Vi{-\b) 

which for 1 < j < k < I form a subset of the whole KP hierarchy in the bilinear form. 



(A21) 



(A22) 



The leading term as e, — > in ( A20 ) corresponds to (j, fc, I) = (1, 2, 3) in ( A22 ). In this case eq. ( |A22| ) 
gives (the bilinear form of) the KP equation itself: 



(£>£ - 4D tl D t3 +3Dl) T -T = 0. 



(A23) 



This example shows once again that the discrete hierarchy has a more transparent structure than the 
continuous one. The continuum limit brings artificial complications. 
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